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Abstract 

This  thesis  investigates  a  cross  range  smear  phenomenon  seen  in  far  field,  fre¬ 
quency  domain  Xpatch®  calculations.  The  phenomenon  is  very  subtle,  manifesting 
itself  in  2-D  ISAR  images  as  a  low-level  scatterer  response  smeared  in  the  cross 
range  direction.  This  cross  range  smear  occurs  only  using  complex  target  models 
with  certain  characteristics.  It  is  also  a  mathematical  construct,  not  occurring  in 
physical  SAR  systems.  Using  a  carefully  constructed  scattering  target  set,  Xpatch®- 
generated  ISAR  images  are  used  to  characterize  cross  range  smear  in  terms  of  its 
input  parameters.  The  characterization  is  done  as  a  DOE-based  polynomial  approx¬ 
imation  to  the  observed  smear  levels.  Frequency  extent  and  bandwidth  have  the 
highest  effect  on  cross  range  smear,  consistently  increasing  smear  with  parameter 
value.  Ray  density  is  slightly  less  important,  having  primarily  squared  and  second 
order  influence.  The  choice  of  diffraction  and  first  bounce  algorithm  has  very  little 
effect  on  cross  range  smear.  In  addition,  the  performance  of  a  proposed  smear  reduc¬ 
tion  technique  is  analyzed  against  Xpatch®-generated  ISAR  images.  The  algorithm 
generally  reduces  smear,  but  the  smear  reduction  magnitude  is  not  a  linear  function 
of  smear  value. 
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CROSS  RANGE  SMEAR  CHARACTERIZATION  IN  XPATCH 


ISAR  IMAGES 


I.  Introduction 

1.1  Background 

There  have  been  great  advances  in  the  field  of  stealth  design  due  to  the  advent 
of  Computational  Electromagnetics  (CEM).  Incorporating  electromagnetic  stealth 
into  modern  weapon  system  design  enables  great  gains  in  lethality  and  survivability. 
These  stealth  aircraft  penetrate  deep  into  highly  defended  areas  unseen  and  deliver 
ordinance  onto  previously  unreachable  targets. 

CEM  tools  are  commonly  used  to  simulate  electromagnetic  scattering  from 
targets,  a  key  step  in  designing  stealth  aircraft.  This  scattered  field  is  defined  as  the 
reflections  off  a  target  from  incident  electromagnetic  waves-mathematically  stated 
as  the  total  field  in  the  presence  of  the  target  minus  the  field  with  no  target  present. 
The  incident  waves  are  characterized  as  from  a  radar,  friendly  or  adversarial.  By 
simulating  the  electromagnetic  scattering  from  the  target,  radar  signatures  can  be 
developed.  These  signatures  are  used  to  analyze  the  target  for  radar  vulnerabilities 
to  be  exploited  or  protected  against,  i.e.,  the  application  of  stealth  design. 

Inverse  Synthetic  Aperture  Radar  (ISAR)  images  are  often  used  as  a  two- 
dimensional  (2-D)  target  signature.  ISAR  processing  provides  a  scattering  center 
map  within  2-D  space:  range  and  cross  range.  A  variation  of  spotlight  Synthetic 
Aperture  Radar  (SAR),  ISAR  processes  multiple  target  aspects  from  target  rota¬ 
tion  against  a  stationary  radar  platform  into  a  2-D  scattering  strength  image.  The 
ISAR  image  allows  radar  vulnerability  analysis  and  automatic  target  recognition  and 
classification. 
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Xpatch®  is  a  popular  CEM  tool  used  to  simulate  and  visualize  electromagnetic 
scattering  from  target  models,  developed  and  distributed  by  Science  Applications 
International  Corporation  (SAIC)  and  sponsored  by  government  organizations  such 
as  the  Air  Force  Research  Laboratory.  Xpatch®  can  generate  a  variety  of  target 
signatures  from  this  scattering  data,  such  as  ISAR  images.  Xpatch®  is  used  in 
many  US  Government  and  industry  applications  [1]. 

1.2  Problem 

In  certain  ISAR  image  classes  generated  with  Xpatch®,  a  low-level  smear  is 
apparent  in  the  cross  range  direction.  An  example  of  this  smearing  is  illustrated  in 
Figs.  1.1  and  1.2.  Figure  1.1  shows  an  ISAR  target  image  when  constrained  to  a  low 
dynamic  range,  where  only  strong  scattering  contributions  are  shown.  Figure  1.2 
shows  the  same  target  image  constrained  to  a  high  dynamic  range,  including  weaker 
responses  masked  in  the  low  dynamic  range  image.  While  more  overall  scatterer 
response  is  evident  in  the  high  dynamic  range  image,  as  expected,  an  additional 
smearing  effect  in  the  cross  range  (vertical)  direction  is  evident.  This  cross  range 
smear  is  the  subject  of  research. 

The  smearing  is  very  subtle,  manifesting  itself  in  the  ISAR  image  as  an  elon¬ 
gated  scatterer  response.  It  occurs  at  a  very  low-level;  the  smearing  can  only  be 
observed  in  the  SAR  images  when  viewed  with  a  high  dynamic  range,  showing  re¬ 
sponses  many  orders  of  magnitude  less  than  the  overall  peak  response  in  the  image. 

Increasingly  sensitive  systems  require  images  with  a  high  dynamic  range.  These 
systems  can  distinguish  low  level  responses  from  background  noise.  Hence,  low  level 
simulation  errors  become  increasingly  important  problems  when  accurately  simulat¬ 
ing  and  modeling  potential  targets.  Cross  range  smear,  an  anomalous  artifact  in 
simulated  data,  therefore  represents  an  obstacle  to  be  overcome  in  order  to  support 
these  newer  systems. 
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Figure  1.1:  An  ISAR  image  illustrating  the  2-D 
scatterer  response  from  a  target,  given  a  30  dB  dy¬ 
namic  range  and  incident  ray  density  of  10  rays/ A.  No 
cross  range  smear  is  observed.  From  [2j. 


range 


Figure  1.2:  An  ISAR  image  illustrating  the  2-D 
scatterer  response  from  a  target,  given  a  55  dB  dy¬ 
namic  range  and  incident  ray  density  of  10  rays/ A. 
Cross  range  smear  is  seen  as  elongated  responses  in 
the  vertical  direction.  From  [2]. 
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Xpatch®  uses  two  entirely  different,  user  selectable,  methods  when  forming 
ISAR  images:  frequency  or  time  domain  calculations.  While  the  results  should  be 
substantially  equivalent,  cross  range  smearing  is  only  evident  in  ISAR  images  formed 
through  frequency  domain  calculations.  ISAR  images  generated  with  time  domain 
processing1  do  not  show  the  same  smearing  effect. 

Open  literature  documents  the  cross  range  smear  present  in  Xpatch®  ISAR 
images.  Bhalla  and  Ling  describe  the  effect  in  a  1997  IEEE  article,  where  they  at¬ 
tribute  the  smear  to  angular  data  scintillation  noise  [2],  However,  a  robust  cross 
range  smear  characterization  and  complementary  rationale  for  its  presence  is  hereto¬ 
fore  lacking.  This  research  performs  a  characterization  of  cross  range  smear  present 
in  Xpatch®  frequency  domain  ISAR  images.  This  characterization  is  described  in 
terms  of  target  geometry,  electromagnetic  (EM)  physics  simulation  technique,  and 
ISAR  image  generation  parameters. 

1.3  Previously  Proposed  Solution 

According  to  the  Shooting  and  Bouncing  Rays  (SBR)  technique  used  to  cal¬ 
culate  scattered  field  data  (the  precursor  to  ISAR  images)  in  Xpatch®,  a  dense  grid 
of  ray  tubes  is  launched  at  the  scattering  target  from  each  ISAR  imaging  position. 
Bhalla  and  Ling  call  the  initial  location  and  ray  tube  positioning  the  “shooting 
grid  [2].”  Each  shooting  grid  position  is  located  at  a  unique  angle  relative  to  the 
target,  which  translates  to  a  unique  spoke  in  the  spatial  frequency  domain  annulus2. 
At  each  new  position,  no  matter  how  small  the  position  change,  the  grid  moves  to 
a  corresponding  angular  location.  This  repositioning  causes  the  grid-launched  ray 
tubes  to  initially  impact  the  target  at  different  points  during  each  ray  shoot.  The 
change  in  initial  reflection  locations  is  inevitable  and  happens  regardless  of  shooting 
grid  density  or  angular  resolution.  The  target  complexity  coupled  with  the  different 

1Xpatch®  time  domain  processing  is  performed  by  calculated  scattered  field  data  from  a  single, 
central  aspect  and  extrapolating  the  results  over  all  other  aspects  in  the  SAR  aperture. 

2The  concept  of  the  spatial  frequency  domain  annulus  is  described  in  Chapter  III. 
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initial  reflection  points  results  in  a  vastly  different  scattered  field,  even  for  small  an¬ 
gular  changes  in  the  shooting  grid.  Bhalla  and  Ling  identified  these  sharp  scattered 
field  variations  (a  function  of  incident  ray  angle)  as  angular  scintillation  noise. 

Bhalla  and  Ling  propose  two  solutions  to  minimize  the  presence  of  cross  range 
smear.  The  first,  a  spatial  filtering  algorithm,  attenuates  high-frequency  scintillation 
noise.  The  second  uses  angular-extrapolation  to  fill  in  and  smooth  gaps  between  the 
discretely  sampled  SBR-calculated  fields. 

1.4  Research  Overview  &  Chapter  Layout 

The  purpose  of  this  research  is  to  characterize  cross  range  smear  in  Xpatch® 
ISAR  images.  This  characterization  has  two  facets.  First,  cross  range  smear  is 
characterized  in  terms  of  Xpatch®  input  parameters  and  target  geometries.  Second, 
a  solution  proposed  by  Bhalla  and  Ling,  described  in  Sec.  1.3,  is  evaluated  using  a 
number  of  Xpatch®-generated  ISAR  images. 

Chapter  II  presents  an  analytic  development  of  key  underlying  concepts.  The 
SBR  technique,  used  by  Xpatch®  to  calculate  scattered  fields,  is  described  with  its 
various  input  parameters.  A  mathematical  development  of  ISAR  follows,  describing 
the  concepts  behind  the  transformation  of  scattered  field  data  into  a  2-D  radar 
image.  The  last  section  of  this  chapter  details  the  algorithm  used  to  determine  the 
amount  of  cross  range  smear  contained  in  an  ISAR  image,  necessary  to  perform  the 
subsequent  smear  characterization  and  reduction  evaluation. 

Chapter  III  follows  the  analytical  development  with  a  methodology  description. 
The  geometries  and  parameter  values  used  to  generate  Xpatch®  ISAR  images  are 
described.  The  techniques  used  to  characterize  smear  in  terms  its  parameters  are 
also  detailed,  as  are  the  smear  reduction  application  and  subsequent  performance 
analysis. 

The  results  of  both  research  facets,  smear  characterization  and  smear  reduc¬ 
tion  analysis,  are  split  into  two  chapters.  Chapter  IV  presents  the  results  of  smear 
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characterization  in  terms  of  each  target  geometry,  as  well  as  a  general  discussion  of 
the  aggregate  results.  Chapter  V  presents  the  performance  of  the  smear  reduction 
technique,  again  in  terms  of  individual  geometries  (but  multiple  ISAR  images,  given 
varying  input  parameters)  and  as  an  aggregate  discussion. 

Finally,  Chapter  VI  presents  final  conclusions.  A  very  general  discussion  of  the 
results  is  given,  along  with  direction  for  further  related  research. 

1.5  Assumptions 

Cross  range  smear  is  a  concern  due  to  the  assumption  that  it  is  a  simulation 
artifact,  a  result  of  the  scattering  algorithm.  This  assumption  implies  cross  range 
smearing  would  not  occur  in  ISAR  images  generated  from  measured  data  given  an 
appropriately  low  noise  floor. 

This  primary  assumption  also  forms  the  basis  for  this  research.  If  smearing 
were  a  physical  effect  present  in  ISAR  images  of  measured  data,  Xpatch®  would  cor¬ 
rectly  model  the  important  scattering  physics  involved.  The  scattered  data  would 
be  correctly  simulated  (at  least  in  terms  of  the  phenomenon  that  produces  the  cross 
range  smear)  and  there  would  be  no  need  to  characterize  the  smearing  as  a  mathe¬ 
matical  anomaly.  Two  primary  reasons  justify  the  assumption  of  cross  range  smear 
existing  as  only  a  simulation  artifact. 

First,  cross  range  smear  generated  in  Xpatch®  utilizing  time  domain  calcula¬ 
tions  does  not  exist.  Since  the  time  domain  scattering  and  ISAR  processing  imple¬ 
mentation  differs  significantly  from  the  frequency  domain  implementation  and  cross 
range  smear  does  not  exist  in  images  generated  using  the  time  domain  processing, 
one  (or  both)  of  the  processing  methods  is  incorrect. 

Second,  ISAR  images  generated  from  measured  data  do  not  exhibit  cross  range 
smear.  Since  no  simulation  is  involved  when  the  scattered  fields  are  directly  mea¬ 
sured,  simulation-specific  problems  are  removed  as  an  error  source.  Assuming  noise, 
present  in  all  measured  signals,  does  not  mask  cross  range  smear  in  measurement- 


6 


based  ISAR  images,  the  assumption  that  cross  range  smears  in  Xpatch®  ISAR  im¬ 
ages  are  a  simulation  artifact  is  justified. 


7 


II.  Analytical  Development 


2. 1  Introduction 

This  chapter  presents  a  description  of  the  techniques  used  both  by  the  Xpatch® 
program  and  the  tools  used  to  analyze  its  output  for  cross  range  smear.  The  Shooting 
and  Bouncing  Rays  (SBR)  technique,  used  by  Xpatch®  to  calculate  the  scattered 
fields  from  a  target  illuminated  by  incident  electromagnetic  (EM)  energy,  is  the  first 
covered.  The  data  calculated  via  SBR  is  then  used  as  inputs  into  Inverse  Synthetic 
Aperture  Radar  (1SAR)  processing,  described  next.  The  ISAR  processing  generates 
a  two-dimensional  (2-D)  image  of  the  target  initially  radiated  by  the  SBR  technique; 
cross  range  smear,  the  subject  of  analysis,  is  sometimes  visible  within  these  ISAR 
images.  Interactions  between  the  SBR  technique  and  ISAR  processing  are  discussed 
after  both  are  presented  individually.  Finally,  a  development  of  the  algorithm  used 
to  analyze  ISAR  images  for  smear  is  presented.  This  algorithm  is  used  to  help  assess 
which  input  factors  influence  the  amount  of  smear  present  in  an  ISAR  image. 

2.2  Shooting  And  Bouncing  Rays 

SBR  is  the  scattering  simulation  method  used  in  Xpatch®,  the  program  that 
generates  ISAR  images  that  can  display  cross  range  smear.  The  SBR  technique  is 
used  to  generate  the  expected  scattered  return  from  a  given  illumination,  such  as  by 
a  simulated  radar.  SBR  is  derived  from  the  combination  of  two  well  known  high- 
frequency  computational  techniques:  Geometrical  Optics  (GO)  and  Physical  Optics 
(PO). 

Balduaf  et  al  described  the  SBR  technique  in  1991  when  they  presented  a  paper 
comparing  SBR  to  measured  data  for  several  standard  scattering  targets  [3].  SBR 
compared  favorably  to  the  measured  data,  especially  with  multiple-bounce  targets 
at  high  frequencies. 


In  SBR,  a  densely  packed  grid  of  GO-type  ray  tubes  is  launched  towards  a 
CAD-generated  target  to  calculate  the  scattered  field  when  illuminated  from  that 
point.  Each  ray  tube  behaves  according  to  the  rules  of  GO  until  its  last  reflection. 
Physical  target  geometries  that  affect  ray  divergence,  convergence,  and  polarization 
are  all  taken  into  account.  At  each  ray  tube  reflection  point,  a  PO  surface  current 
is  calculated.  The  computed  surface  currents  are  then  used  (by  reradiating  those 
currents)  to  find  the  resulting  scattered  field. 

For  the  case  of  cavity  analysis,  a  slightly  different  approach  is  used  [4],  Each 
ray  tube  that  enters  the  cavity  is  reflected  as  before.  However,  the  last  reflection 
point  is  defined  as  the  cavity  opening.  Aperture  integration  is  used  to  calculate  the 
“cavity  equivalent”  of  surface  currents.  The  scattered  field  due  to  these  equivalent 
currents  is  then  calculated. 

Each  ray  tube  individually  contributes  to  a  unique  surface  or  equivalent  cur¬ 
rent.  Therefore,  the  entire  scattered  field  must  be  updated  for  each  ray  tube  launched 
at  the  target. 

In  the  following  sections,  Xpatch®  parameters  are  individually  discussed. 

2.2.1  Ray  Density.  As  mentioned  in  Sec.  2.2,  the  individual  GO  ray  tubes 
propagated  towards  the  target  are  packed  together  into  a  densely  packed  grid.  This 
grid,  referred  to  as  the  shooting  grid,  is  positioned  at  each  incident  radiation  direc¬ 
tion.  The  shooting  grid  defines  the  initial  position  of  ray  tubes  to  be  propagated 
towards  the  target  at  each  shooting  grid  position.  The  grid  span  is  large  enough 
to  cover  the  entire  projected  target  area,  considering  the  target  orientation  and  il¬ 
lumination  direction.  The  ray  tube  density  within  the  shooting  grid  is  called  the 
grid  density.  The  higher  the  grid  density,  in  rays  per  wavelength  (A),  the  more 
ray  tubes  are  propagated  towards  the  target  from  each  position.  Each  ray  tube  is 
spaced  regularly  over  the  shooting  grid.  An  example  shooting  grid  layout  is  shown 
in  Fig.  2.1. 
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Figure  2.1:  An  illustration  of  the  shooting  grid  from  which  ray  tubes 
propagate  towards  a  target  in  SBR  scattering  calculations.  The  shooting 
grid  defines  the  initial  position  of  the  ray  tubes,  and  in  this  case  has  a 
length  of  2A,  height  of  A,  and  grid  density  of  5  rays  per  A. 

2.2.2  First  Bounce  Algorithm.  The  first  bounce  algorithm  parameter  mod¬ 
ifies  the  default  behavior  of  Xpatch® ’s  use  of  the  SBR  technique  when  calculating 
scattered  fields.  The  Xpatch®  first  bounce  parameter  has  three  options:  SBR,  PO, 
and  Z-Buffer.  The  first  two  options,  SBR  and  PO,  are  discussed.  The  Z-buffer  op¬ 
tion  is  not  considered  as  it  provides,  according  to  the  Xpatch®  documentation,  the 
least  accurate  result. 

When  SBR  is  used  as  the  first  bounce  algorithm,  incident  rays  are  launched 
at  the  target  from  the  direction  of  incidence.  This  behavior  is  consistent  with  the 
behavior  of  SBR  at  subsequent  points  of  reflection.  If  target  model  surfaces  are 
illuminated  by  a  ray  tube,  a  reflection  is  generated  at  that  surface.  However,  if 
PO  is  used  as  the  first  bounce  algorithm,  each  target  model  facet  or  patch  is  split 
into  subsections.  Visibility  from  these  subsections  to  rays  from  the  direction  of 
incidence  is  determined  and  PO  currents  are  calculated.  PO  provides  a  more  accurate 
first  bounce  surface  current  model  than  SBR  with  an  accompanying  increase  in 
computational  time. 
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2.2.3  Maximum.  Bounces.  The  maximum  bounces  parameter  is  used  in  the 
Xpatch®  SBR  implementation  to  determine  the  number  of  reflections  that  each  ray 
tube  is  allowed  before  truncation.  This  parameter  is  applicable  to  ray  tubes  that  are 
reflected  many  times  before  reradiating  off  the  target  into  free  space.  Limiting  the 
number  of  reflections  can  decrease  computation  time  at  the  cost  of  scattered  field 
accuracy. 

2.2.4  Contributing  Bounces.  Xpatch® ’s  SBR  implementation  allows  for 

selection  of  which  bounces  contribute  to  the  resulting  scattered  field.  If  the  option 
for  all  bounces  to  contribute  to  the  scattered  field  is  selected,  every  reflection  off  a 
surface  generates  a  surface  current  and,  therefore,  a  scattered  field.  These  fields  are 
then  added  into  the  total  scattered  field  result.  This  option  is  the  most  accurate. 
Alternatively,  only  the  first  and  last  bounces  can  be  set  to  contribute  to  the  total 
scattered  field.  This  option  reduces  computation  time,  at  the  expense  of  accuracy. 

2.2.5  Higher  Order  Bounces.  Xpatch®  can  selectively  choose  which  fields 
contribute  to  the  surface  currents  at  each  SBR  bounce  point  on  the  illuminated 
target.  The  surface  currents  are  then  used  to  calculate  the  resulting  total  scattered 
field.  This  parameter  is  applicable  when  all  SBR  reflections  are  set  to  contribute  to 
the  total  scattered  field,  as  in  Sec.  2.2.4.  Either  both  the  scattered  and  incident  field 
or  only  the  incident  field  at  each  bounce  point  can  be  selected  to  contribute.  When 
both  the  incident  and  scattered  field  contribute,  both  fields  are  used  to  determine 
the  surface  currents  at  each  bounce  point.  If  only  the  scattered  field  is  selected,  only 
the  scattered  field  is  used  to  determine  the  surface  currents.  The  selection  of  both 
incident  plus  scattered  fields  produces  the  most  accurate  result. 

2.2.6  Divergence  Factor.  Xpatch®  can  use  the  divergence  factor,  a  function 
of  surface  curvature,  to  allow  for  ray  tube  spreading  from  bounces  on  curved  surfaces. 
SBR  ray  divergence  affects  the  cross  sectional  area  of  each  ray  tube,  and  therefore 
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the  resulting  scattered  field  from  each  bounce.  This  behavior  is  identical  to  how 
GO  ray  divergence  is  handled.  This  parameter  does  not  apply  to  facetized  target 
models  since  surface  curvature  is  not  defined.  However,  for  target  models  defined  in 
formats  such  as  the  Initial  Graphics  Exchange  Specification  (IGES),  ray  divergence 
is  a  factor. 

2.2.7  Diffraction.  Xpatch®  can  use  a  diffraction  edge  file  to  separately 
calculate  scattered  EM  fields  due  to  specified  diffractive  edges.  Edge  diffraction 
is  calculated  using  Mitzner’s  incremental  edge  diffraction  coefficient  method,  and 
is  included  only  when  this  option  is  specified  and  the  edge  is  unobstructed.  The 
diffractive  scattered  field  is  coherently  added  to  the  scattered  field  calculated  using 
SBR. 

2.3  Inverse  Synthetic  Aperture  Radar 

ISAR,  or  inverse  SAR,  is  the  process  used  by  Xpatch®  to  generate  2-D  radar 
images  of  a  given  scattering  target.  By  using  data  calculated  by  the  SBR  technique 
in  Sec.  2.2,  Xpatch®  uses  ISAR  processing  to  generate  a  2-D  radar  image  of  the 
target  scene.  These  images  can  display  the  smearing  phenomena  that  is  the  subject 
of  research. 

ISAR  is  a  variation  of  spotlight  mode  SAR,  a  well  established  and  very  effective 
technique  that  provides  high-resolution  imagery  of  a  spatial  area  by  leveraging  radar 
return  information  from  multiple  target  aspects.  Beamsteering  is  typically  used 
to  continually  focus  radar  energy  on  a  single  target  area,  thereby  generating  the 
multiple  aspect  information.  ISAR,  as  simulated  in  Xpatch®,  reverses  the  spotlight 
mode  scenario  by  rotating  the  target  while  using  a  stationary  radar.  Target  rotation 
generates  the  multiple  aspects  required  by  ISAR  processing. 

Munson  established  the  basis  for  current  spotlight  mode  processing  techniques  [5]. 
Given  the  similarities  between  spotlight  mode  SAR  and  ISAR,  these  techniques  ap- 
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ply  to  both  processing  types.  He  adapted  techniques  and  theorems  commonly  used  in 
computer-aided  tomography,  better  known  as  CAT  scans,  to  simplify  spotlight  mode 
processing.  In  particular,  he  described  the  processing  in  terms  of  the  projection-slice 
theorem. 

Jakowatz  characterized  the  projection-slice  theorem  by  stating  that  the  one¬ 
dimensional  (1-D)  Fourier  transform  of  a  linear  projection  of  the  reflectivity  density 
(the  target  area  imaged)  is  equal  to  the  2-D  Fourier  transform  of  the  reflectivity 
density  [6].  Mathematically,  the  projection-slice  theorem  is  written  as 

G(U cos#,  U sin#)  =  Pg(U),  (2.1) 


where  G  is  the  2-D  Fourier  transform  of  the  reflectivity  density,  Pq  is  the  Fourier 
transform  on  the  linear  projection,  and  9  is  the  angle  at  which  the  linear  projection  of 
the  reflectivity  density  is  taken.  All  of  these  functions  are  in  the  U  domain,  defined 
by  a  d-rotated  coordinate  system  in  the  spatial  frequency  domain.  This  equation 
allows  the  projection  function  (derived  from  the  SAR  pulse  return)  to  be  related 
to  the  underlying  reflectivity  density  (the  resulting  SAR  image).  The  1-D  Fourier 
transform  of  the  linear  projection  function,  pg,  from  the  spatial  (u)  domain  into  the 
spatial  frequency  (U)  domain,  is 


Pe(U)=  /  pe(u)e-juUdu 


(2,2) 


and  the  2-D  Fourier  transform  of  the  reflectivity  density  g(x,y)  is 


G(U  cos  6,  U  sin  9)  =  G(X,  Y) 


t*  oo  r  oo 


g(x,y)e  rtxX+vY) dxdy . 


'  — oo  J  — OO 


(2,3) 


Jakowatz  further  described  how  a  Linear  Frequency  Modulated  (LFM)  radar 
pulse  can  be  used  in  conjunction  with  the  projection-slice  theorem  to  gather  data 
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for  spotlight  mode  images.  The  LFM  pulse  is  given  by 


s(t) 


cos  (c u0t  +  at2)  n 


(2.4) 


where  ujq  is  the  pulse  center  frequency,  a  is  the  ramping  constant,  and  fl  is  the 
rectangular  envelope  function  that  evaluates  to  unity  between  ±rc/2  (where  rc  is  the 
total  signal  time  length)  and  zero  elsewhere. 

The  incident  LFM  pulse  is  scattered  by  the  target  reflectivity  density  g  with  a 
spatial  convolution  of  the  emitted  LFM  pulse  with  linear  projections  of  the  reflectiv¬ 
ity  density  orthogonal  to  the  direction  of  signal  propagation.  The  linear  projections, 
used  to  represent  the  simultaneous  sampling  of  reflectivity  density  along  a  line  de¬ 
fined  by  the  signal’s  planar  wave  front,  are  given  in  the  spatial  ( u )  domain  as 

/OO 

g[x(u,v),y(u,v)\dv.  (2.5) 

-OO 


where  u  and  v  are  defined  by  x  =  ucosO  and  y  =  usind.  This  relationship  is  shown 
in  Fig.  2.2.  The  LFM  pulse  return  is  then  written  as 

n(t)  =  A  yu‘  pe(u)s  ( t  -  2-A±A)  (2.6) 

where  A  is  an  attenuation  constant  due  to  distance,  R  is  the  distance  to  the  reflectiv¬ 
ity  density  center,  c  is  the  speed  of  light,  the  LFM  pulse  s(t)  is  given  by  Eqn.  (2.4), 
and  the  integration  limits  represent  the  cross  range  extent  of  the  linear  projection. 

The  LFM  pulse  return,  Eqn.  (2.6),  is  deramped  after  reception.  Deramping 
is  a  process  where  the  incoming  signal,  in  this  case  the  radar  return,  is  quadrature 
demodulated  (to  I  and  Q  channels)  and  low  pass  filtered.  The  deramped  channels 
are  then  recombined  into  one  signal.  The  first  step  to  deramping  is  to  mix  the 
received  signal  with  inphase  and  quadrature  signals  delayed  to  match  the  two-way 
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Figure  2.2:  An  illustration  of  the  linear  projection  pg(u )  on 
a  ^-rotated  ( u ,  v)  coordinate  system  over  a  reflectivity  density 
g(x,y). 

propagation  time  of  the  received  signal.  The  inphase  mixing  term  is  given  by 


C/(t)  =  cos 


,  2  R\  (  2  R\ 

w°l  t-T)+a{t-T) 


and  the  quadrature  mixing  term  is 


cQ(t)  =  sin 


,  2  R\  (  2  R 

WO  P - —  Q  / - 

c  )  \  c 


(2.7) 


(2.8) 


Alixing  each  these  terms  with  the  received  signal  in  Eqn.  (2.6)  and  utilizing  the 
trigonometric  identities 


1 

cos  A  cos  B  =  -  [cos(5  —  A)  +  cos(A  +  B )]  (2.9) 

cos  A  sin  B  =  -  [sin(£>  —  A)  +  sin(A  +  B)\  (2-10) 
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results  in  processed  signals  given  by 


re, lit)  = 


f*ll  i 


Pe(u) 


'  -u  1 


42-o(t-f-^)+a(i-f)2+a(t-™)2; 


+e' 


n  l^~2(i?  +  M)/cS)  du 

Tr 


(2.11) 


and 


?e,Q(t)  = 


A, 

2 


ru  i 


-3< 


r  —  u± 


D,[c0^+a(t-f)2-a(t-m); 


+e' 


n(^-2(^  +  n)/cM^, 

r 


(2.12) 


where  3?  and  S  take  the  real  or  imaginary  parts  of  the  integral,  respectively,  to 
extract  the  appropriate  sine  or  cosine  terms.  The  tilde  (~)  on  the  received  signals 
denotes  the  received  signals  have  been  changed  by  processing. 

The  next  deramping  step  is  low  pass  filtering.  Low  pass  filtering  these  two 
mixed  signals  removes  the  second  exponential  term,  the  sum  term  in  Eqns.  (2.9)  and 
(2.10).  Slightly  rearranging  the  first  exponential  term  allows  writing  the  low  pass 
filtered  signals  as 


h,,(t)  = 


ru  i 


Po(u 


’  ~U\ 


)e^'ha(^)2+(^)(^0+2Q(*“^))]n  ( t~2(R  +  u)/c 

V 


du 


(2.13) 


and 


reM  =  ^ 


ru  i 


Pe(u 


'  —U  i 


)elHw)2+(^)(-°+2K*-w))]n  ( t~2(i?  +  n)/c 

\  7~r 


du 
(2.14) 
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These  I  and  Q  terms  are  then  combined  to  form 


h{t)  = 


Pe(u ) 


,4-Q(v)2+(f)(-o+2 


n 


t  —  2(R  +  u)/c 


Tr 


du 


(2.15) 


Ignoring  the  first  term  in  the  exponential,  — ck(^)2,  and  attributing  it  to  skew  allows 
rewriting  Eqn.  (2.15)  as 


h{t) 


p0(M)eJ'[w(^°+2Q(t-¥))]  n 


-  2 (R  +  u)/c^j 


(2.16) 


This  equation  can  now  be  recognized  as  the  Fourier  transform  of  po,  and  can  be 
rewritten  as 


ro{t) 


JP«{-C  (^o  +  2a 
^MU) 


(2.17) 

(2.18) 


where  t  has  been  time  limited  by  II.  The  effect  of  this  time  limiting  is  discussed 
below.  By  Eq.  2.17,  the  spatial  frequency  U  domain  maps  to  the  time  domain  as 


U 


2 

c 


+  2ck 


(2.19) 


where  t  is  again  time  limited  by  II. 

According  to  the  projection  slice  theorem  shown  in  Eqn.  (2.1),  Pq  represents  a 
linear  trace  through  G  (the  2-D  Fourier  transform  of  the  reflectivity  density).  Since, 
the  U  domain  is  a  d-rotated  version  of  the  spatial  frequency  domain,  the  deramped 
return  can  be  viewed  per  the  projection  slice  theorem  as 

fe(t)  =  ^G(X,Y)  (2.20) 

=  —G(Ucos9, 1/sind).  (2.21) 
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Each  processed  return,  at  a  unique  9,  traces  a  spoke  into  G(U cos  9,U sin#). 
Multiple  returns  at  different  9  values  form  a  disc  sector  G.  Time  limiting  the  return 
rg(t)  constrains  the  disc  to  an  annulus,  where  annulus  width  is  proportional  to  LFM 
pulse  bandwidth  and  the  offset  is  proportional  to  the  center  pulse  frequency.  This 
relationship  is  illustrated  in  Fig.  2.3. 

From  G(U  cos  9,  [/sin  9),  the  2-D  Fourier  transform  of  the  reflectivity  den¬ 
sity,  a  SAR  image  can  be  reconstructed  by  a  2-D  inverse  Fourier  transform.  Since 
G(U  cos  9 ,  U  sin  9)  is  sampled  on  the  polar  grid  naturally  defined  by  the  spokes  shown 
in  Fig.  2.3,  a  polar  to  rectangular  interpolation  is  usually  performed.  This  inter¬ 
polation  then  allows  Fast  Fourier  Transform  (FFT)  use.  The  2-D  inverse  Fourier 
transform  projects  the  data  from  the  spatial  frequency  ( U )  domain  to  the  image  (u) 
domain,  defined  in  range  and  cross  range. 


2.3.1  Bandwidth.  Bandwidth  defines  the  total  frequency  change  the  tar¬ 
get  is  illuminated  through.  Usually  measured  in  Gigahertz  (GHz),  the  bandwidth 
directly  determines  the  ISAR  image  range  resolution.  This  range  resolution  <5range  is 
found  from 

Grange  =  (2.22) 

where  c  is  the  speed  of  light,  B  is  the  propagated  signal  bandwidth,  and  <5range  is  the 
ISAR  image  range  resolution. 

Since  bandwidth  is  inversely  proportional  to  resolution,  a  more  detailed  image 
can  generally  be  generated  using  a  higher  bandwidth.  Frequency  extent  determines 
the  width  of  the  spatial  frequency  arc  shown  in  Fig.  2.3,  and  both  frequency  increment 
and  bandwidth  together  define  the  number  of  frequencies  that  are  considered  for  each 
target  aspect. 

The  bandwidth  maps  to  the  spatial  frequency  annulus  shown  in  Fig.  2.3  as 


,  AnB 

annulus  width  =  - , 

c 


(2.23) 
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Figure  2.3:  Spotlight  data  gathered  in  spatial  frequency  space. 
The  angle  of  the  data  in  spatial  frequency  is  the  same  as  the  angle 
at  which  the  data  is  physically  collected.  Annulus  width  is  propor¬ 
tional  to  the  radar  pulse  bandwidth,  and  offset  is  proportional  to 
pulse  center  frequency. 
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where  B  is  the  bandwidth  of  the  signal  in  Hz ,  and  c  is  the  speed  of  light.  This 
mapping  corresponds  to  Eqn.  2.19,  where  the  second  term  determines  the  annulus 
width  and  can  therefore  be  written  as 


annulus  width 
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(2.26) 

(2.27) 

(2.28) 


where  t  is  evaluated  at  its  fl  extents  and  the  bandwidth  B  is  related  to  the  pulse 
ramping  constant  o;  by  B  =  (q>/ti)tc  [7]. 


2.3.2  Frequency  Center.  The  frequency  center  defines  the  midpoint,  in 
frequency,  of  the  signal  illuminating  the  target.  This  frequency  center  also  determines 
how  far  the  spatial  frequency  arc,  shown  in  Fig.  2.3,  is  displaced  from  the  origin. 
The  higher  the  frequency  center,  the  shorter  the  wavelength  propagated  towards  the 
target. 

The  frequency  center  maps  to  the  spatial  frequency  annulus  shown  in  Fig.  2.3 
as 

47T 

annulus  midpoint  =  — ,  (2.29) 

^0 
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where  Ao  is  the  mean  (or  center)  pulse  wavelength.  This  mapping  corresponds  to 
Eqn.  2.19,  where  lvq  is  the  frequency  center  and 


annulus  midpoint 
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(2.30) 

(2.31) 

(2.32) 

(2.33) 


The  first  term  in  Eqn.  2.19  thus  determines  the  midpoint  of  the  annulus,  and  relates 
to  the  frequency  center. 


2.3.3  Frequency  Increment.  Since  returned  signal  frequencies  are  discretely 
sampled,  the  frequency  increment  defines  the  separation  between  consecutive  signal 
frequencies.  The  frequency  increment  also  defines  the  total  range  extent  of  the 
resulting  ISAR  image  as 
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(2.34) 

(2.35) 


where  <5/  is  the  frequency  increment,  B  is  the  bandwidth  in  Hz,  and  Arange  is  the 
total  ISAR  image  range  extent.  Frequency  increment  and  bandwidth  together  define 
the  number  of  frequencies  the  received  signal  is  evaluated  for.  Each  frequency  thus 
has  a  returned  magnitude  and  phase  associated  with  it. 


2.3.4  Angular  Span.  The  angular  span  defines  the  angle  over  which  the 
radar  illuminates  the  target.  The  angular  span  is  equal  to  the  angular  span  of  the 
spatial  frequency  arc,  Fig.  2.3,  extends.  Assuming  a  narrowband  signal,  the  angular 
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span  affects  the  resulting  ISAR  image  resolution  as 


5 


crossrange 
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(2.36) 


where  A0  is  the  mean  or  center  propagated  signal  wavelength,  Ae  is  the  total  angle 
span  the  signals  are  propagated  from,  and  <5Crossrange  is  the  resulting  cross  range 
resolution.  Angular  increment  and  angular  span  together  determine  the  number  of 
aspects  from  which  a  signal  must  be  propagated  towards  a  target. 


2.3.5  Angular  Increment.  Angular  increment  defines  the  separation,  in 
angle,  between  consecutive  propagation  locations.  The  angular  increment  defines 
regular  intervals  inside  the  total  angular  span  defined  by  Sec.  2.3.4  from  which  sig¬ 
nal  frequencies  (set  by  the  frequency  span  and  center  of  Secs.  2.3.1  and  2.3.2)  are 
propagated.  The  resulting  ISAR  image’s  cross  range  extent  Acrossrange  is 
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where  6$  is  the  angular  increment,  A g  is  the  total  angular  span,  and  Ao  is  the  average 
or  center  wavelength  of  the  propagated  signal.  Given  an  angular  center,  angular 
increment  and  angular  span  determine  the  aspects  from  which  signals  are  propagated 
towards  the  target. 


(2.37) 

(2.38) 


2-4  SBR  and  ISAR  Interaction 

Xpatch®  uses  SBR  to  calculate  the  scattered  fields  from  a  given  set  of  aspects 
and  illumination  frequencies.  The  scattered  fields  are  used  by  the  ISAR  process¬ 
ing  algorithm  to  generate  the  appropriate  ISAR  image.  While  the  SBR  and  ISAR 
processing  algorithms  are  distinct,  the  parameters  of  each  can  affect  the  behavior 
of  the  other.  This  section  describes  the  interaction  between  SBR  and  ISAR  when 
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Figure  2.4:  The  process  that  Xpatch®  uses  to  generate  an  ISAR  image 
starts  with  the  definition  of  target  geometry,  SBR-specific  parameters,  and 
ISAR-specific  parameters.  This  data  feeds  into  the  SBR  process,  which 
generates  the  scattered  fields  at  the  specified  aspects  and  frequencies.  The 
scattered  fields  are  then  processed  into  an  ISAR  image. 


Xpatch®  generates  an  ISAR  image  from  the  parameters  specified  in  Secs.  2.2  and 
2.3.  An  illustration  of  SBR  and  ISAR  process  and  parameter  interaction  is  shown 
in  Fig.  2.4.  While  both  SBR  and  ISAR  parameters  influence  the  final  scattered 
fields  and  resulting  ISAR  images,  only  SBR  parameters  change  the  SBR  algorithm’s 
inherent  behavior. 

Parameters  specific  to  SBR,  covered  in  Sec.  2.3,  propagate  directly  through 
the  ISAR  process.  Each  parameter  influences  the  calculated  scattered  field.  Since 
the  scattered  field  is  input  to  the  ISAR  image  generation  process,  these  parameters 
affect  the  final  ISAR  image. 

Parameters  specific  to  ISAR,  covered  in  Sec.  2.3,  form  the  requirements  dictat¬ 
ing  the  aspect  and  frequency  set  from  which  the  SBR  technique  is  applied.  Although 
they  affect  the  application  of  the  SBR  algorithm,  the  SBR  algorithm  itself  is  not  al¬ 
tered. 


2-4-1  Ray  Density.  Ray  density,  described  in  Sec.  2.2.1,  affects  the  scat¬ 
tered  field  calculation  by  reflecting  at  specific  points  over  the  target  geometry.  In 
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true  EM  scattering,  the  incident  wave  is  usually  described  as  a  continuous  plane 
wave  over  the  entire  target1.  However,  the  inherent  discretization  of  ray  tubes  in 
the  SBR  technique  is  only  an  approximation  of  an  incident  plane  wave.  The  inci¬ 
dent  ray  tubes  impinge  discrete  points  on  the  target  geometry.  The  ray  tubes  then 
successively  reflect  onto  another  set  of  discrete  points  on  the  geometry  until  they 
are  reradiated  into  free  space.  A  different  incident  aspect  angle  produces  a  different 
set  of  initial  reflection  points.  This  different  initial  reflection  point  sets  can  result 
in  dramatically  different  reflection  angles,  greatly  impacting  the  final  scattered  field 
and  the  resulting  ISAR  image. 

The  scattered  field  maps  into  the  spatial  frequency  domain  for  ISAR  processing 
and  is  then  Fourier  transformed,  producing  the  final  image.  Sharp  variances  with 
respect  to  aspect  in  the  spatial  frequency  domain  are  synonymous  with  frequency 
content  throughout  the  Fourier  transformed  image  domain.  Since  the  polar  to  rect¬ 
angular  interpolation  performed  before  the  Fourier  transform  substantially  preserves 
discontinuities,  any  discontinuities  in  the  spatial  frequency  angular  direction  result 
in  image  content  spread  throughout  the  cross  range  extent  of  the  resulting  ISAR 
image.  This  content  is  in  the  same  range  bin  as  the  discontinuity  that  caused  it. 
Aspect  dependent  discontinuities  have  a  greater  effect  on  cross  range  image  content 
than  range  content  due  to  mapping  constant  aspect  frequency  sweeps  along  annulus 
spokes  (as  in  Fig  2.3). 

Aspect-based  discontinuities  due  to  the  inherent  shooting  grid  discretization 
are  most  pronounced  for  geometries  formed  of  shapes  changing  significantly  with 
respect  to  angle.  For  example,  a  cavity  produces  a  dramatically  different  scattered 
held  with  a  small  incident  aspect  change.  This  result  occurs  due  to  the  large  num¬ 
ber  of  reflections  the  small  change  in  initial  reflection  points  propagates  through. 
Similarly,  a  target  geometry  with  sharp  edges  produces  a  dramatically  different  set 

lrThis  description  is  of  course  also  an  approximation,  relying  on  the  assumption  that  the  target 
is  completely  enclosed  within  a  quiet  zone  where  the  incident  spherical  wave  is  locally  planar 
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of  initial  reflection  points  with  a  small  change  in  incident  aspect,  attributed  to  the 
incident  ray  reflecting  on  opposite  sides  of  the  edge  as  the  aspect  changes.  The  ray 
tube  then  continues  until  finally  reflected  off  geometry.  The  final  reflection  point 
may  be  significantly  displaced  from  the  last  reflection  point  of  the  same  ray  tube  at 
the  previous  aspect.  A  flat  plate,  on  the  other  hand,  does  not  change  the  initial  in¬ 
cident  reflection  points  much  based  on  a  small  change  in  incident  aspect.  Therefore, 
the  resulting  scattered  field  changes  are  minimal. 

2-4-2  Angular  Increment.  Angular  increment,  described  in  Sec.  2.3.5,  is  a 
parameter  that  determines  the  final  ISAR  image’s  cross  range  dimension.  Decreas¬ 
ing  the  aspect  increment  results  in  smaller  aspect  changes  between  each  subsequent 
SBR  ray  shoot.  Aspect  increment  does  not  affect  the  SBR  process  itself,  but  it  does 
affect  the  incident  aspect  locations  where  SBR  is  applied.  Therefore,  any  change 
in  aspect  increment  alone  does  not  affect  the  scattered  held  calculation;  each  cal¬ 
culation  is  independently  applied.  However,  aspect  increment  variation  does  have 
a  significant  affect  on  the  final  ISAR  image.  The  smaller  the  aspect  increment,  the 
more  spokes  are  present  in  the  spatial  frequency  domain.  This  increased  spatial  fre¬ 
quency  resolution  results  a  higher  rectangularly  resampled  resolution  and  leads  to  a 
correspondingly  higher  cross  range  image  extent. 

2-4-3  Angular  Span.  Angular  span,  described  in  Sec.  2.3.4  and  used  with 
the  angular  increment  parameter,  determines  the  discrete  set  of  angles  where  the 
incident  SBR  rays  are  propagated.  While  the  angular  span  controls  the  final  ISAR 
image  cross  range  resolution,  it  does  nothing  to  inherently  change  the  SBR  algorithm 
behavior.  The  SBR  algorithm  is  simply  applied  over  an  aspect  defined  in  part  by 
angular  span. 

2-4-4  Frequency  Center.  Frequency  center,  described  in  Sec.  2.3.2  and  used 
in  combination  with  frequency  increment  and  bandwidth  to  determine  the  discrete 
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frequency  set  the  SBR  calculated  scattered  fields  are  evaluated.  Once  the  points  of 
reflection  and  path  lengths  are  calculated  for  a  given  aspect,  the  scattered  fields  are 
evaluated  in  closed  form  for  each  frequency.  Altering  the  frequencies  does  not  change 
the  ray  reflection  angle  or  magnitude,  only  its  resultant  path  length  phase.  Since 
all  rays  are  separately  propagated  and  the  scattered  fields  from  each  propagated  ray 
are  summed  in  SBR,  the  SBR  calculated  total  scattered  held  varies  as  a  function 
of  frequency.  This  variation  naturally  leads  to  variation  in  the  final  ISAR  image 
content. 

2-4-5  Frequency  Increment.  Frequency  increment,  described  in  Sec.  2.3.3 
and  used  in  combination  with  frequency  center  and  bandwidth  determines  the  dis¬ 
crete  frequency  set  SBR  calculated  scattered  fields  are  evaluated.  Once  the  points 
of  reflection  and  path  lengths  are  calculated  for  a  given  aspect,  the  scattered  fields 
are  evaluated  in  closed  form  for  each  frequency.  Therefore  only  one  ray  shoot  is 
performed  at  each  aspect,  regardless  of  the  number  of  frequencies  at  that  aspect. 

More  frequencies  are  evaluated  at  each  aspect  as  the  frequency  increment  is 
reduced.  This  reduction  causes  Xpatch®  to  evaluate  more  scattered  held  values, 
mapped  along  each  spoke  in  the  spatial  frequency  domain.  In  addition  to  increasing 
the  image  extent  of  the  hnal  ISAR  image,  better  rectangular  to  polar  interpolation 
can  be  performed  prior  to  the  Fourier  transform  that  creates  the  ISAR  image. 

2-4-6  Bandwidth.  Bandwidth,  described  in  Sec.  2.3.1  and  used  in  com¬ 
bination  with  frequency  center  and  frequency  increment  to  determine  the  discrete 
frequency  set  SBR  calculated  scattered  helds  are  evaluated.  Since  an  entire  ray 
shoot  is  only  performed  once  at  each  aspect,  increasing  the  bandwidth  does  not  add 
ray  shoots  or  change  the  performance  of  the  SBR  algorithm  itself.  This  increase 
will,  however,  affect  the  hnal  ISAR  image  resolution  by  providing  a  wider  spatial 
frequency  domain  annulus. 
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2.5  Smear  Assessment  Algorithm 

In  order  to  measure  the  input  parameter  and  geometry  effects  on  cross  range 
smear,  a  quantitative  method  assessing  the  cross  range  smear  extent  is  required. 
This  assessment  algorithm  mathematically  generates  a  numerical  value  for  the  smear 
present  in  an  ISAR  image.  This  numerical  smear  value  feeds  into  the  analysis  of 
parameters  affecting  smear  in  Xpatch®  ISAR  images. 

Since  ISAR  images  are  essentially  2-D  magnitude  plots,  worst  case  smear  man¬ 
ifests  itself  as  constant  real- valued  magnitude  across  the  entire  cross  range  extent. 
In  other  words,  the  cross  range  slice  under  examination  consists  of  one  constant 
value  across  the  entire  cross  range.  Conversely,  minimum  smear  is  a  purely  real 
delta  function  within  the  cross  range  slice.  The  delta  functions  represent  the  scat¬ 
tering  centers  of  the  imaged  reflectivity  density.  Further  discussion  on  maximum 
and  minimum  smear  is  given  as  the  smear  assessment  algorithm  is  evaluated  for 
each  case. 

First,  a  single  cross  range  slice  g(x,yo)  is  taken  from  the  image.  Due  to  ISAR 
image  pixelization,  the  slice  is  defined  discretely  along  x.  The  slice  is  then  trans¬ 
formed  into  its  frequency  domain  using  a  one  dimensional  Fourier  transform  and 
designated  as  G(X,y0).  The  result  is  not  in  a  true  spatial  frequency  domain  since 
a  2-D  Fourier  transform  is  required  in  the  image  formation  process,  so  the  result  is 
denoted  G  instead  of  G, 


G(X,y0)  =  X{g(x,y0)}. 


(2.39) 


The  absolute  value  of  G( X,  y0)  is  taken  to  extract  its  magnitude, 


\G(X,y0)\  =  \X{g(x,y0)}\ , 


(2.40) 
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and  then  differentiated  over  X  as 

^{\G(X,y0)\}  =  -^{\^{g(x,y0)}\}.  (2.41) 

This  differentiation  is  discrete  because  of  the  discrete  nature  of  X ,  so  the  derivative 
is  simply  the  difference  between  consecutive  values  of  \G(X,  y0)|.  Since  cross  range 
smear  is  a  result  of  angle-based  discontinuities  in  the  spatial  frequency  data,  this 
step  provides  an  approximate  measure  of  discontinuity. 

An  absolute  value  is  again  applied  to  force  all  decreases  in  \G(X,y0)\,  and 
hence  all  negative  values  of  ^\G(X,y0)\,  positive.  This  step  is  applied  as 

T{|G(J>f,!/o)l}  =  ■  (2.42) 

Now  that  all  changes  in  \G(X,  t/0)|  contribute  positively,  each  value  of  the 
discrete  function  ^\G(X,  y0)|  is  summed  to  produce  the  single  smear  value.  This 
smear  value  can  take  any  value  in  [0,  oo) ,  where  minimum  smear  is  when  the  func¬ 
tion  \G{X,  2/0)  |  is  constant  and  maximum  smear  is  approached  when  \G(X,  y0)|  has 
large  changes/discontinuities.  The  smear  value  in  each  cross  range  slice  is  therefore 
constructed  as 

smear (7/0)  =  sum  j  J-  ||G(X,j/0)||  j  (2.43) 


28 


=  sum 


(2.44) 


The  smear  value  associated  with  the  entire  ISAR  image  is  the  maximum  of  all 
cross  range  slice  smear  values,  or 

=  max  | sum  j  ^  {\F {g(x,  j/0)}|}  |  J  •  (2.46) 

2.5.1  Analysis  of  Maximum  Smear.  Extensive  smear  in  a  magnitude  ISAR 
image  is  a  constant  real- valued  slice  extending  the  entire  cross  range  length.  Further, 
maximum  smear  is  present  when  the  smear  is  constant  over  an  ISAR  image  infinite 
in  cross  range  extent.  The  constant  cross  range  slice  is  then  denoted  as 

g(x,y0)  =  C.  (2.47) 

Its  Fourier  transform  is  a  purely  real  delta  function, 

F{C}  =  C6(X).  (2.48) 

Taking  the  absolute  value  has  no  effect  since  everything  is  already  positive,  so 

\T{C}\  =  \C5(X)\  (2.49) 

=  C5(X).  (2.50) 


(2.45) 
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The  derivative  of  this  function,  implemented  as  the  difference  between  consecutive 
discrete  X  values,  is 

|Ctf(X)|  (2-51) 

— CS(X )  (2.52) 

ax 

oo  X  =  0— 

<  -oo  X  =  0+  (2.53) 

0  else 

where  the  infinite  values  are  located  at  the  intervals  (— <5_y,0)  and  (0,  Sx)  given  the 
X  resolution  Sx-  Another  absolute  value  ensures  all  the  derivatives  are  positive,  in 
order  to  constructively  contribute  to  the  sum.  The  above  result  then  becomes 

X  =  0- 

X  =  0+  (2.54) 

else. 

The  summation  over  all  X ,  defined  as  the  smear  in  this  cross  range  slice  with  maxi¬ 
mum  smear,  approaches  infinity. 

As  a  practical  example,  the  Matlab®-implemented  cross  range  smear  algorithm 
is  tested  on  a  unity  constant  over  a  2500  sample  cross  range  slice.  This  input  function 
is  shown  in  Fig.  2.5.  In  this  case,  the  algorithm  assesses  a  13440  smear  value. 

2.5.2  Analysis  of  Minimum  Smear  With  Single  Point  Scatterer.  A  reflec¬ 
tivity  density  can  generally  be  decomposed  into  a  discrete  and  independent  point 
scatterer  set,  called  scattering  centers  [8].  These  scattering  centers  are  ideally  repre¬ 
sented  in  an  ISAR  image  as  delta  functions  in  range  and  cross  range.  Therefore,  an 
ISAR  image  without  smear  is  a  discrete  2-D  delta  function  set.  Each  delta  function 
represents  a  scattering  center  of  the  reflectivity  density  being  imaged. 
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Figure  2.5:  A  cross  range  smear  slice  with  maxi¬ 
mum  smear,  using  2500  cross  range  pixels,  generates  a 
13440  smear  value. 

A  single  scattering  center  illustrates  a  best  case  scenario,  an  ISAR  image  with¬ 
out  smear.  Since  this  single  point  scatterer  is  shown  as  a  delta  function  in  the  ISAR 
image  domain  and  the  ISAR  image  presents  no  phase  information,  the  scattering 
center  is  given  as  a  delta  function  located  at  xn  scaled  by  a  purely  real  constant  C. 
As  the  image  is  analyzed  over  each  cross  range  slice,  assumed  infinite  in  extent,  the 
scattering  center  presents  itself  as 


g(x,y0)  =  CS(x  -  xn).  (2.55) 

Its  Fourier  transform  is  a  purely  real  constant  appropriately  frequency  shifted, 

F{C8{x)}  =  Ce~jXnX.  (2.56) 
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where  X  is  finite  in  extent  and  discretely  sampled.  Taking  the  magnitude  of  this 
value  strips  off  the  complex  exponential, 

\F{C8{x)}\  =  \Ce~iXnX\  (2.57) 

=  C.  (2.58) 


The  derivative  of  this  constant  is  zero  for  all  discrete  values  of  X, 


d 

dx 


?{C5{x)}  | 


dx 

0. 


(2.59) 

(2.60) 


This  result  is  summed  over  all  X  to  produce  a  smear  value  for  a  single  point  scatterer 
of  identically  zero. 

As  a  practical  example,  the  Matlab®  implemented  cross  range  smear  algorithm 
was  tested  on  a  single  delta  functions  of  unity  strength  displaced  70%  down  a  cross 
range  slice  with  2500  samples.  This  input  function  is  shown  in  Fig.  2.6.  In  this  case, 
the  algorithm  assesses  a  smear  value  of  4.19e-13. 

2.5.3  Analysis  of  Minimum  Smear  With  Multiple  Point  Scatterers.  The 
result  for  minimum  smear  can  be  generalized  for  a  discrete  set  of  independent  point 
scatterers  along  the  cross  range  slice.  This  set  represents  an  accurate  point  scatterer 
decomposition  of  a  reflectivity  density  imaged  into  a  magnitude- only  ISAR  plot 
as  a  set  of  real  delta  functions  arbitrarily  placed  along  the  cross  range  slice.  The 
corresponding  reflectivity  density  is  given  by 

g{x,y0)  =  ^2,Cn5{x  -  xn),  (2.61) 

n 

where  xn  is  the  spatial  location  and  Cn  is  the  positive  magnitude  of  delta  function 
n.  The  cross  range  extent  is  assumed  infinite.  Its  Fourier  transform  is  the  sum  of 
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Figure  2.6:  A  single  delta  function  displaced  20% 
from  the  center  of  a  2500  pixel  cross  range  slice  gener¬ 
ates  a  smear  value  of  4.19e-13. 


frequency-shifted  constants, 


T 


Y,CJ(x-xn)  >  =  Cne~jx 


(2.62) 


where  X  is  discretely  sampled  and  finite  in  extent.  Taking  the  magnitude  of  this 
value  removes  the  phase  of  the  result.  However,  since  the  complex  exponentials 
constructively  and  destructively  interfere,  the  magnitude  will  smoothly  vary  within 
bounds  set  by  [0,E nCn\.  This  step  is  shown  as 


T 


(2.63) 


The  variation  between  consecutive  points  of  Eqn.  (2.63)  will  depend  on  the  number 
of  point  scatterers  n  and  their  position  in  the  cross  range  slice  xn.  The  greater  n  is, 
the  sharper  variation  between  consecutive  points  in  Eqn.  (2.63)  can  be. 

The  derivative  of  Eqn.  2.63  is  the  difference  between  values  at  consecutive 
points  along  X.  The  value  at  each  point  along  X  of  this  derivative  will  depend  on 
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n  and  xn.  However,  since  the  differentiated  function  is  the  sum  of  sinusoids  and 
therefore  without  large  discontinuities,  the  difference  between  consecutive  points  in 
Eqn.  (2.63)  is  kept  small  (but  non-zero). 

Taking  an  absolute  value  ensures  that  all  discontinuities  are  represented  by 
positive  numbers.  This  result  is  then  summed  over  the  X  extent  to  cumulatively 
measure  all  discontinuities.  Since  each  difference  is  relatively  small,  their  sum  will  be 
much  less  than  an  otherwise  identical  maximum  smear  case.  In  the  maximum  smear 
case,  dramatic  discontinuities  contribute  very  large  differences  to  the  summation. 

As  a  practical  example,  the  Matlab®  implemented  cross  range  smear  algorithm 
was  tested  on  three  delta  functions  of  varying  strength  placed  arbitrarily  along  a  cross 
range  slice  with  2500  samples.  This  input  function  is  shown  in  Fig.  2.7.  In  this  case, 
the  algorithm  assesses  a  smear  value  of  2513.  This  number  is  much  less  than  the 
13440  computed  for  the  maximum  smear  case  of  Fig.  2.5. 

Another  example  is  shown  in  Fig.  2.8.  Here,  five  delta  functions  of  varying 
strength  are  placed  arbitrarily  along  the  same  cross  range  slice.  This  function  has 
a  smear  value  of  3019,  which  is  still  significantly  less  than  the  13440  computed  for 
maximum  smear. 
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Figure  2.7:  Three  delta  functions  placed  along  a 
cross  range  slice  of  2500  samples  have  a  smear  value  of 
2513. 


Figure  2.8:  Five  delta  functions  placed  along  a  cross 
range  slice  of  2500  samples  have  a  smear  value  of  3019. 


III.  Methodology 


3. 1  Introduction 

This  chapter  presents  the  methodologies  used  to  perform  two  tasks:  character¬ 
ize  cross  range  smear  and  analyze  the  smear  reduction  method  presented  by  Bhalla 
and  Ling  [2].  The  smear  characterization  is  given  in  terms  of  the  geometries  and  in¬ 
put  parameters  used  by  Xpatch®  when  generating  Inverse  Synthetic  Aperture  Radar 
(ISAR)  images.  The  smear  reduction  implementation  and  analysis  methodology  is 
then  described,  building  upon  the  smear  characterization  work. 

Cross  range  smear,  present  in  Xpatch®  ISAR  images  and  the  subject  of  char¬ 
acterization,  is  a  product  of  scattered  fields  generated  through  Xpatch® ’s  Shooting 
and  Bouncing  Rays  (SBR)  implementation.  Chapter  II  describes  the  SBR  technique 
and  ISAR  image  generating  process.  The  methodology  characterizing  cross  range 
smear  is  separated  into  the  following  tasks. 

1.  Identify  parameters  potentially  affecting  cross  range  smear. 

2.  Analytically  describe  how  parameters  affect  cross  range  smear. 

3.  Define  scattering  targets  with  geometries  inducing  a  variety  of  scattering  mech¬ 
anisms. 

4.  Define  a  set  parameter  values  and  target  geometries,  used  to  generate  Xpatch® 
ISAR  images,  that  allows  important  cross  range  smear  factors  to  be  extracted 
and  identified. 

5.  Generate  Xpatch®  ISAR  images  with  chosen  targets  and  variable  parameters. 

6.  Analyze  the  factors  (parameters  and  geometries)  affecting  cross  range  smear 
in  the  resulting  Xpatch®  ISAR  images. 

The  parameters  affecting  SBR  and  subsequent  ISAR  image  generation  were  previ¬ 
ously  defined  and  individually  discussed  in  Secs.  2.2  and  2.3.  Section  2.4  further 
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discusses  these  parameters  in  terms  of  their  interaction  with  the  SBR/ISAR  process 
and  potential  for  affecting  cross  range  smear. 

The  scattering  targets  and  parameter  extents  discussed  in  this  chapter  are  used 
to  generate  Xpatch®  ISAR  images.  Section  3.2  describes  the  various  targets  used 
to  generate  each  ISAR  image  set.  For  each  scattering  target,  an  ISAR  image  set 
is  generated  by  varying  previously  identified  parameters.  Section  3.3  identifies  the 
specific  parameter  values  and  ranges.  A  Design  of  Experiments  (DOE)  approach 
is  used  to  construct  an  ISAR  image  set  using  specific  parameter  values  facilitating 
an  empirical  cross  range  smear  analysis  as  a  function  of  those  input  parameters. 
Section  3.4  describes  the  DOE  analysis  approach  while  Chap.  IV  presents  the  results 
of  this  analysis. 

This  chapter  also  discusses  the  implementation  of  the  method  presented  by 
Bhalla  and  Ling  [2]  to  reduce  cross  range  smear.  Section  3.5  describes  how  the 
smear  reduction  is  implemented  and  how  the  subsequent  smear  analysis  is  performed. 
Chapter  V  presents  the  results  of  this  analysis.  In  addition  to  the  smear  reduction 
technique  as  proposed  by  Bhalla  and  Ling,  a  slightly  modified  version  is  presented. 

Xpatch®  version  4.7.16  is  used  for  all  scattered  held  computation  and  ISAR 
image  generation. 

3.2  Target  Geometries 

This  section  identifies  target  models  used  by  Xpatch®  along  with  parameters 
defined  in  Sec.  3.3  when  calculating  the  scattered  fields.  ISAR  images  are  generated 
from  these  scattered  holds  and  then  analyzed  for  cross  range  smear  presence.  These 
targets  are  Computer-Aided  Design  (CAD)  models,  each  with  geometric  features 
affecting  electromagnetic  (EM)  scattering.  This  scattering  forms  the  foundation  for 
the  ISAR  imaging  process,  therefore  careful  geometry  selection  is  extremely  impor¬ 
tant.  The  targets  chosen  are  also  of  varying  complexity,  generally  divided  into  two 
groups:  canonical  and  moderately  complex  targets. 
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Table  3.1:  The  target  geometries  exercise  a  variety  of 
scattering  mechanisms  in  Xpatch®  while  generating  scat¬ 
tered  field  data  and  the  resulting  ISAR  image. 


Geometry  Name 

Category 

Type 

Flat  Rectangular  Plate 
Ridged  Plate 

Slicy 

Specular 

Edge 

Multiple 

Canonical 

Moderately  Complex 
Moderately  Complex 

Canonical  shapes  provide  one  dominant  scattering  mechanism,  such  as  specular 
reflection,  edge  diffraction,  or  multiple  (2-3)  bounce  reflections.  Moderately  complex 
targets  possess  several  scattering  mechanisms  affecting  the  resultant  scattered  field, 
such  as  cavity-type  multiple  bounce  reflections  or  additional  diffraction  components. 
A  table  of  the  target  geometries  used  is  given  in  Table  3.1. 

Each  geometry  file  is  described  in  one  of  two  CAD  formats:  IGES  or  facet.  An 
IGES  file  exactly  describes  curved  surfaces  while  a  facet  file  describes  flat  triangular 
facets  that,  at  best,  approximate  a  curved  surface.  Each  respective  section  specifies 
the  CAD  file  format  used  for  each  geometry. 

Diffraction  due  to  geometric  edges  is  defined  in  Xpatch®  by  a  separate  (.edge) 
file.  This  file  identifies  diffraction-inducing  straight  lines  of  the  CAD  geometry. 
Since  diffraction  is  not  considered  in  purely  SBR  calculated  scattering,  the  diffractive 
edges’  scattering  contribution  is  only  added  to  the  primary  SBR  scattering  when  the 
parameter  to  include  diffraction  effects,  Sec.  3. 3. 1.3,  is  enabled. 

Other  common  scattering  mechanisms,  such  as  creeping  and  traveling  waves, 
are  not  modeled  due  to  the  nature  of  the  SBR  high  frequency  assumptions  and 
application. 

All  target  geometries  are  imaged  with  a  waterline  ISAR  aperture.  The  SAR 
aperture,  or  the  collection  of  imaging  aspects  defined  by  its  angular  span  and  an¬ 
gular  increment,  is  thus  set  along  zero  degrees  elevation  and  sweeps  out  an  angular 
span  centered  about  zero  degrees  azimuth  (centered  on  the  x-axis)  as  illustrated  by 
Fig.  3.1. 
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Figure  3.1:  An  illustration  of  a  waterline  ISAR  aperture,  at  zero  degrees 
elevation  and  centered  about  zero  degrees  azimuth  (the  x  axis).  The  ISAR 
aperture  in  this  example  extends  for  a  total  angular  span  of  twenty  degrees. 
All  ISAR  images  generated  for  geometries  in  Sec.  3.2  use  this  general  orien¬ 
tation. 
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Figure  3.2:  The  flat  rectangular  plate  geometry  induces 
primarily  single-bounce  specular  scattering. 

3.2.1  Flat  Rectangular  Plate.  The  flat  rectangular  plate  geometry  shown  in 
Fig.  3.2  is  a  100  inch  by  100  inch  plate  normal  to  the  x  axis,  with  a  uniform  thickness 
of  1  inch.  As  a  primarily  specular  geometry,  the  effects  due  to  SBR  application  are 
exclusively  single  bounce  reflections.  Diffraction  along  all  12  edges,  4  on  each  side 
plus  4  due  to  the  width,  is  also  present.  However,  diffraction  is  secondary  to  the 
primary  specular  scattering. 

The  flat  rectangular  plate  serves  as  a  canonical  shape  that  examines  the  SBR 
and  ISAR  interactions  due  to  primarily  specular  scattering.  It  is  specified  as  an 
IGES  file. 

3.2.2  Ridged  Plate.  The  ridged  plate  geometry  shown  in  Fig.  3.3  is  a  flat 
rectangular  plate  with  triangular  ridges  over  the  entire  +x  surface,  extending  in  the 
z  direction.  The  1  inch  thick  base  plate  is  100  inches  by  100  inches  in  extent.  Each 
ridge  is  an  additional  10  inches  from  base  to  tip  and  is  10  inches  along  each  base  for 


40 


Figure  3.3:  The  ridged  plate  geometry  is  a  flat  plate  with 
triangular  ridges  facing  in  the  positive  x  direction.  It  is 
mostly  a  single  and  double  bounce  scattering  geometry  with 
edge  diffraction  along  the  edges  of  the  base  plate  and  on  the 
top  of  each  ridge.  The  ridges  are  meant  to  excite  a  large  an¬ 
gular  reflection  change  and,  therefore,  a  large  change  to  the 
scattered  field  based  on  a  small  incident  angular  increment. 
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Figure  3.4:  The  slicy  geometry  is  a  composite  of  several 
scattering  geometries  including  specular  surfaces,  singly 
curved  surfaces,  dihedrals,  trihedrals,  and  cavities.  These 
geometries  induce  specular  and  multiple  bounce  reflections, 
as  well  as  diffraction  along  exterior  edges. 

a  total  of  10  ridges  over  the  entire  +x  surface.  Diffraction  is  due  to  the  edges  along 
the  plate,  as  well  as  the  top  of  each  ridge.  Scattering  on  the  ridges  is  primarily  due 
to  single  and  double  bounce  reflections.  The  ridges  induce  dramatic  changes  in  the 
ray  reflections  and  scattered  field  due  to  a  small  change  in  incident  angle. 

The  ridged  plate  is  a  medium  complex  shape  that  examines  the  SBR  and 
ISAR  interactions  due  to  primarily  low-order  reflections.  Those  reflections  vary 
considerably  with  incident  angular  direction.  It  is  specified  as  an  IGES  hie. 

3.2.3  Slicy.  The  slicy  geometry  shown  in  Fig.  3.4  includes  many  different 
scattering  mechanisms  such  as  cavities,  cylinders,  dihedrals,  trihedrals,  and  singly 
curved  surfaces.  Its  extents  are  defined  by  a  bounding  box  that  is  96  inches  along  the 
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x  direction,  108  inches  in  the  y  direction,  and  66  inches  in  the  z  direction.  Diffraction 
is  from  all  exterior  edges  with  a  wedge  (interior)  angle  of  90  degrees. 

The  slicy  geometry  is  a  medium  complex  shape  that  examines  SBR  and  ISAR 
interactions.  This  geometry  induces  many  different  scattering  mechanisms.  It  is 
specified  as  a  facet  hie,  so  ray  divergence  is  ignored  and  all  curved  surfaces  are 
facetized. 

3.3  Parameters 

Xpatch®  generates  scattered  held  data  and  the  resulting  ISAR  image  using  a 
set  of  input  parameters  along  with  a  specihed  target  geometry.  This  section  discusses 
the  parameter  values  used  to  generate  ISAR  images  that  are  analyzed  for  cross  range 
smear. 

Each  parameter  is  either  varied  or  held  constant.  Those  that  are  varied,  iden¬ 
tified  and  discussed  in  Chap.  II,  potentially  affect  cross  range  smear  in  the  output 
ISAR  images.  The  value  ranges  for  each  of  these  variable  parameters  are  discussed 
and  justified  in  Sec.  3.3.1.  Xpatch®  needs  additional  parameters  to  produce  ISAR 
images,  but  these  do  not  necessarily  influence  cross  range  smear.  These  parameters 
are  held  constant  for  all  ISAR  images  and  are  discussed  and  justified  in  Sec.  3.3.2. 

3.3.1  Variable  Parameters.  This  section  presents  the  parameters  varied 
when  producing  the  ISAR  image  set.  This  image  set  is  then  analyzed  for  cross 
range  smear.  Table  3.2  presents  a  succinct  variable  parameter  list.  Continuous 
parameters  may  take  any  value  between  the  given  maximum  and  minimum  extent 
while  nominal  parameters  cycle  between  explicitly  declared  values.  These  values 
constitute  the  extent  the  parameters  are  allowed  vary  over  in  the  course  of  cross 
range  smear  analysis. 

3. 3. 1.1  Ray  Density.  Ray  density,  analytically  discussed  in  Sec.  2.2.1, 
is  a  continuously  variable  parameter.  It  is  constrained  to  vary  between  5  and  20 
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Tabic  3.2:  Variable  Parameters 


Parameter 

Values 

Units 

Range  Type 

Ray  Density 

5-20 

Rays/A 

Continuous 

First  Bounce  Algorithm 

PO,SBR 

N/A 

Nominal 

Diffraction 

yes,  no 

N/A 

Nominal 

Bandwidth 

1-10 

GHz 

Continuous 

Frequency  Center 

10-15 

GHz 

Continuous 

Frequency  Increment 

0.005-0.01 

GHz 

Continuous 

Angular  Span 

5-20 

degrees 

Continuous 

Angular  Increment 

0.01-0.1 

degrees 

Continuous 

rays/A.  The  Xpatch®  default  is  10  rays/A,  so  a  factor  of  two  is  used  to  define  the 
maximum  and  minimum  ray  density  values. 

3.3. 1.2  First  Bounce  Algorithm.  The  first  bounce  algorithm,  analyt¬ 
ically  discussed  in  Sec  2.2.2,  is  a  nominally  variable  parameter  set  to  either  PO  or 
SBR.  A  third  option,  Z-buffer,  exists,  however  it  is  infrequently  used  and  computa¬ 
tionally  comparable  to  SBR.  SBR  is  the  recommended  Xpatch®  default,  while  PO 
provides  additional  accuracy  at  the  cost  of  computational  efficiency.  Therefore,  both 
SBR  and  PO  are  possible  values  for  this  parameter. 

3.3. 1.3  Diffraction.  Diffraction  is  a  nominally  variable  parameter 
used  by  Xpatch®  to  add  a  diffraction  effect  to  the  the  total  (specular)  scattered  field 
otherwise  calculated  by  the  SBR  technique.  EM  diffraction  is  separately  computed 
using  an  additional  geometry  hie  defining  edge  presence.  The  diffraction  component 
is  then  coherently  added  to  the  fields  calculated  via  the  SBR  technique.  Since  EM 
diffraction  can  have  a  strong  effect  on  the  total  scattered  field,  it  is  an  important 
parameter  to  examine.  The  inclusion  of  diffraction  is  specified  by  a  yes  or  no  value. 


3.3. 1-4  Bandwidth.  Bandwidth,  analytically  discussed  in  Sec.  2.3.1, 
is  a  continuously  variable  parameter.  It  is  constrained  to  vary  between  1  and  10 
GHz. 
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3.3. 1.5  Frequency  Center.  Frequency  center,  analytically  discussed 
in  Sec.  2.3.2,  is  a  continuously  variable  parameter.  It  is  constrained  to  vary  between 
10  and  15  GHz. 

3. 3. 1.6  Frequency  Increment.  Frequency  increment,  analytically  dis¬ 
cussed  in  Sec.  2.3.3,  is  a  continuously  variable  parameter.  It  is  constrained  to  vary 
between  0.005  and  0.01  GHz. 

3. 3. 1.7  Angular  Span.  Angular  span,  analytically  discussed  in  Sec.  2.3.4, 
is  a  continuously  variable  parameter.  It  is  constrained  to  vary  between  5  and  20 
degrees.  Since  angular  span  determines  the  resulting  ISAR  image  cross  range  res¬ 
olution,  this  interval  allows  for  a  sub-unit  cross  range  resolution  fine  enough  to 
distinguish  smaller  geometry  features. 

3.3. 1.8  Angular  Increment.  Angular  increment,  analytically  dis¬ 
cussed  in  Sec.  2.3.5,  is  a  continuously  variable  parameter.  It  is  constrained  to  vary 
between  0.01  and  0.1  degrees.  Since  angular  increment  determines  the  cross  range 
ISAR  image  extent,  this  interval  ensures  that  all  geometries  are  fully  contained  within 
the  ISAR  image  extent. 

3.3.2  Constant  Parameters.  This  section  presents  the  parameters  held 
constant  in  the  generation  of  each  ISAR  image,  along  with  their  rationale  for  use. 
Table  3.3  presents  a  succinct  list  of  these  constant  parameters  along  with  their  ap¬ 
propriate  values.  These  parameters,  when  combined  with  the  variable  parameters 
defined  in  Sec.  3.3.1  and  a  target  geometry  in  Sec.  3.2,  form  the  inputs  Xpatch® 
needs  to  compute  scattered  fields  and  generate  a  resulting  ISAR  image. 

3.3.2. 1  Max  Bounces.  Max  Bounces,  analytically  discussed  in  Sec.  2.2.3, 
is  the  maximum  number  of  ray  bounces  a  SBR  ray  tube  can  reflect  through  before 
being  truncated.  The  Xpatch®  default  is  50.  A  higher  value  only  affects  situations 
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Tabic  3.3:  Constant  Parameters 


Parameter 

Value 

Max  Bounces 
Contributing  Bounces 
Higher  Order  Bounces 
Divergence  Factor 
Computational  Domain 
Signature  Entry  Mode 
Square  Resolution 

50 

all 

Incident  +  Scattered  Field 
yes 

Frequency 

SAR  Image 

no 

where  an  extremely  large  number  of  reflections  occur,  such  as  cavities.  The  scattered 
fields  due  to  such  reflections  can  usually  be  considered  negligible.  Hence,  this  value 
was  left  constant  and  unchanged. 

3. 3. 2. 2  Contributing  Bounces.  The  contributing  bounces  parameter, 
analytically  discussed  in  Sec.  2.2.4,  is  set  to  the  Xpatch®  recommended  value  of 
‘all’.  Due  to  the  computations  required  for  each  reflection,  computing  the  additional 
scattered  field  contributions  from  every  reflection  point  (up  to  the  value  of  Max 
Bounces)  does  not  impose  significant  computational  overhead.  In  addition,  scattered 
fields  are  more  accurately  computed  when  all  bounces  are  taken  into  account. 

3. 3. 2. 3  Higher  Order  Bounces.  Xpatch®  calculates  the  surface  cur¬ 
rents  at  each  reflection  point  using  the  higher  order  bounces  parameter.  When  the 
parameter  is  set  to  ‘Incident  +  Scattered  Field’,  both  the  incident  and  scattered 
held  contributions  at  the  reflection  point,  in  addition  to  the  ray  tube  being  reflected, 
contribute  to  the  surface  current.  This  parameter  setting  provides  the  most  accurate 
surface  current  calculation,  and  the  most  accurate  scattered  field  results.  The  other 
possibility  involves  only  including  the  scattered  field  and  reflected  ray  tube.  How¬ 
ever,  using  both  the  incident  and  scattered  held  with  the  ray  tube  is  recommended. 
This  parameter  is  available  when  Contributing  Bounces,  Sec.  3. 3. 2. 2,  is  set  to  ‘all’. 
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3. 3. 2. 4  Divergence  Factor.  The  divergence  factor  parameter  is  used 
by  Xpatch®  to  change  the  cross  section  of  each  SBR  ray  tube  based  on  the  surface 
curvature  at  the  reflection  point.  It  is  a  standard  GO  method  of  compensating  for 
surface  curvature  to  more  accurately  compute  the  scattered  held.  This  parameter 
is  only  applicable  to  IGES  models,  where  surface  curvature  is  defined.  Facetized 
models  do  not  have  any  surface  curvature.  Therefore,  this  parameter  has  no  effect 
on  those  targets.  Divergence  factor  is  enabled  for  all  cases. 

3. 3. 2. 5  Computational  Domain.  Xpatch®  determines  the  type  of 
scattering  computations  needed  and  the  need  for  additional  parameters  (such  as 
frequencies  and  aspects)  using  the  option  for  selecting  the  computational  domain  in 
tandem  with  the  signature  entry  mode  in  Sec.  3. 3. 2. 6.  Since  cross  range  smear  is 
a  problem  in  frequency  domain  ISAR  images,  the  computational  domain  is  set  to 
Frequency  Domain. 

3. 3. 2. 6  Signature  Entry  Mode.  Xpatch®  determines  whether  the 
scattering  computation  should  be  done  in  the  frequency  or  time  domain,  along  with 
the  primary  scattered  held  hie  format  (.ss  or  .si),  using  the  signature  entry  mode 
parameter.  Since  cross  range  smear  is  generated  through  frequency  domain  com¬ 
putations,  SAR  Image  is  set  as  the  constant  value  for  this  parameter.  Selecting 
this  value  allows  for  the  input  of  other  necessary  items,  such  as  frequency  and  as¬ 
pect  parameters  that  determine  ISAR  image  dimensions,  and  produces  a  complex 
held/spatial  frequency  (.ss)  hie  from  which  an  ISAR  image  is  derived.  Any  other 
value,  such  as  RCS  or  Range  Prohle  does  not  produce  the  necessary  scattering  data 
for  an  ISAR  image. 

3. 3. 2.1  Square  Resolution.  Xpatch®  uses  the  square  resolution  pa¬ 
rameter  to  force  the  range  and  cross  range  resolutions  to  be  equal.  Enabling  this 
parameter  does  not  allow  independent  variation  of  frequency  and  aspect  parameters 
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and  prevents  the  individual  parameter  characterization.  Therefore,  square  resolution 
is  set  to  ‘no’  for  all  cases. 

3-4  Cross  Range  Smear  Characterization 

With  the  scattering  targets  of  Sec.  3.2  and  the  parameters  of  Sec.  3.3,  a  DOE 
approach  identifies  important  factors  contributing  to  cross  range  smear.  The  DOE 
approach  generates  a  set  of  ‘runs’  [9],  each  varying  the  input  parameters  to  produce 
a  separate  ISAR  image.  Each  target  geometry  is  considered  separately,  so  each  set 
of  runs  is  performed  for  each  target  geometry.  Once  Xpatch®  has  produced  the 
requisite  ISAR  image  set,  each  image  is  analyzed  for  the  existence  and  extent  of 
cross  range  smear  per  Sec.  2.5.  These  cross  range  smear  values  are  then  used  to 
extract  the  influence  of  that  geometries’  input  parameters.  This  process  is  repeated 
for  each  geometry,  and  the  results  form  the  characterization  of  cross  range  smear. 

The  DOE  approach  separates  input  parameters  into  two  categories:  contin¬ 
uous  and  nominal.  Continuous  parameters  may  take  any  value  between  a  defined 
maximum  and  minimum  while  nominal  parameters  only  take  discrete  values,  such 
as  ‘yes’  or  ‘no’.  From  the  output  value  analysis,  an  empirical  model  is  built  approxi¬ 
mating  the  SBR  and  ISAR  image  generation  process — a  mathematical  function  that 
produces  a  smear  value  as  a  function  of  its  input  parameters.  Each  input  parame¬ 
ter  influences  the  amount  of  cross  range  smear  with  first,  second,  and  higher-order 
effects.  First-order  effects  are  weighted  parameter  values,  second-order  effects  are  a 
weighted  product  of  two  (possibly  identical)  input  parameters,  and  so  on.  A  first- 
order  main  effects  model  for  three  parameters  with  interaction  terms  can  be  written 
as 


smear  =  /(pi,  P2,  P3)  (3.1) 

=  Co  +  C1P1  +  C2P2  +  C3p3  +  C12PlP2  +  C13P1P3  +  C23p2P3,  (3.2) 
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where  p  is  an  input  parameter  and  c  is  the  corresponding  weight,  or  coefficient. 
Each  input  parameter  is  normalized  to  vary  linearly  between  -1  and  1.  The  overall 
number  of  terms  (first-order,  second-order,  etc.)  in  the  function  equals  the  number 
of  degrees  of  freedom  (DOFs)  in  the  model  and  is  directly  related  to  the  number 
ISAR  images  analyzed.  The  error  between  the  Xpatch®  SBR  and  ISAR  process  and 
the  DOE-prodnced  model  is  also  directly  related  to  the  DOE  design  and  the  number 
of  ISAR  images  analyzed. 

There  are  two  prominent  experiment  design  types  within  the  DOE  framework: 
full  and  fractional  factorial.  A  full  factorial  design  analyzes  2n  runs  by  testing  all 
combinations  of  high  and  low  parameter  extents,  where  n  is  the  number  of  continuous 
and  two  level  (yes-no)  nominal  input  parameters.  This  design  increases  the  number 
of  required  runs  exponentially  with  the  number  of  parameters,  and  therefore  is  not 
a  good  choice  for  the  eight  input  parameters  considered  under  analysis.  Fractional 
factorial  designs  analyze  a  subset  of  full  factorial  runs  and  are  denoted  by  designs 
where  n  is  the  number  of  parameters,  p  is  the  equivalent  number  of  parameters 
that  are  eliminated  by  reducing  the  number  of  runs,  and  R  is  the  resolution.  The 
resolution  gives  an  indication  of  the  order  of  aliased  parameter  interactions,  where 
different  parameter  effects  are  indistinguishable  and  therefore  modeled  together  as 
one  parameter.  Since  a  rR~p  design  only  tests  high  and  low  parameter  values,  at 
best  a  linear  model  can  be  constructed. 

The  number  of  runs  determines  the  DOFs  in  the  resulting  model.  One  degree 
of  freedom  (DOF)  is  used  to  estimate  the  amount  of  error  in  the  model,  another  is 
typically  used  to  estimate  the  intercept  (constant  term)  of  the  model.  The  DOFs  not 
used  to  estimate  the  error  or  intercept  are  used  to  find  the  contributions  of  individual 
(first-order)  or  combined  (second-order)  input  parameters,  the  terms  in  Eqn.  (3.2). 

While  analyzing  Xpatch®  ISAR  images,  where  each  run  can  last  a  few  hours 
or  days,  it  is  important  to  minimize  the  number  of  runs  while  still  keeping  parameter 
aliasing  to  an  acceptable  level.  A  2®r4  fractional  factorial  design  is  used  that  required 
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16  runs  for  8  (continuous  or  two  level  nominal)  input  parameters.  The  resolution 
of  V  means  this  design  allows  no  aliasing,  or  inseparable  combination,  of  first-order 
effects  with  second-order  effects.  First-order  effects  in  a  resolution  V  design  are 
aliased  with  third-order  and  higher  effects.  Therefore,  the  2®r4  model  is  written  as 

smear  =  f(pi,P2,P3,P4,P5,P6,P7,Ps)  (3.3) 

=  Co  +  CiPi  +  C2p-2  +  C3p3  +  C4p4  +  C5p5  +  C6p6  +  C7p7  +  C-sPs 

+C12P1P2  +  c13pip3  +  c23p2P3  +  •  •  •  (3.4) 

where  p\  is  normalized  ray  density,  p2  is  normalized  bandwidth,  etc.  The  first-order 
aliasing  structure  for  this  design  is 


Pi  =>■  Pi  +  PzPiPd  +  P2P3P7  +  P2PPP3  +  P2PbPn 

+P4P5P7  +  P3P5P8  +  P&P7P8  +  higher  —  order  terms.  (3.5) 

P2  =>  P2+  P3P4P5  +  P1P3P7  +  PiPaPs  +  P1P5P6 

+P4Pep7  +  P3PePs  +  PbP7p8  +  higher  —  order  terms.  (3.6) 

P3  =>  P3+  P2P4P7,  +  PlPm  +  PlP2P7  +  P1P5P8 

+P2P6P8  +  P4P7P8  +  PbP&P7  +  higher  —  order  terms.  (3.7) 

PA  =>  P4  +  P2P3P5  +  P1P3P6  +  P1P2P8  +  P1P5P7 

+P2P6P7  +  P3P7P8  +  P5P&P8  +  higher  —  order  terms.  (3.8) 

P5  =>  P5+  P2P3PA  +  P1P2PQ  +  PlPAP7  +  PlP3P8 

+P3P&P7  +  PaPgPs  +  P2P7P8  +  higher  —  order  terms.  (3.9) 

Pe  =>  P&  +  P1P3PA  +  P1P2P5  +  P2PAP7  +  P2P3P8 

+P3P5P7  +  PaPtPPs  +  P1P7P8  +  higher  —  order  terms.  (3.10) 
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P7  =>  P7+  PlP2P3  +  P\p4p5  +  P2P4P6  +  P3P4P8 

+P3P5P6  +  P2P5P8  +  P1P6P8  +  higher  -  order  terms.  (3.11) 
P8  =>  P8+  P1P2P4  +  PlP'.iPb  +  P2P3P6  +  P3P4P7 

+P4P5P6  +  P2P5P7  +  P1P6P7  +  higher  —  order  terms.  (3.12) 

The  2*7 4  design  and  its  16  runs  allows  14  contributions  to  be  estimated.  These 
contributions  are  a  combination  of  first  and  second-order  parameter  effects.  Unfor¬ 
tunately,  not  all  first  and  second-order  contributions  can  be  estimated.  Since  each 
scattering  geometry  is  considered  separately,  a  different  model  is  constructed  for 
each  model.  The  estimated  contributions  are  therefore  different  for  each  scatter¬ 
ing  geometry;  the  contributions  are  changed  as  part  of  the  analysis  to  best  fit  the 
experimental  cross  range  smear  data  gathered  for  each  model. 

The  16  runs  only  test  parameter  values  at  the  outer  extents  of  the  ranges  of 
each  continuous  parameter  with  combinations  of  the  two  level  nominal  parameters. 
Therefore  the  model  derived  from  these  tests  is  at  best  linear.  By  adding  runs  that 
test  values  at  the  center  point  of  the  continuous  parameters  and  for  each  combination 
of  two  level  nominal  parameters,  the  analysis  can  account  for  response  curvature 
between  the  parameter  extents.  Since  there  are  (2)  two  level  nominal  parameters, 
four  additional  runs  are  used  that  test  the  midpoints  of  each  continuous  parameter. 
With  these  additional  center  point  runs,  the  total  of  runs  is  increased  to  20  and  the 
model  can  account  for  quadratic  terms.  The  smear  model  is  then  written  as 

smear  =  f(pi,P2,P3,P4,P5,P6,P7,Pa)  (3-13) 

=  C0  +  C1P1  +  C2P2  +  C-3P3  +  C4p4  +  C5p5  +  C6p6  +  C7p7  +  C8P8 
+C11P1  +  C12P1P2  +  C13P1P3  +  C14P1P4 

+C22P2  +  C23P2P3  +  C24P2P4  +  •  •  •  (3-14) 

The  aliasing  structure  is  unaffected  by  added  center  point  runs. 
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The  JMP  statistical  software  package,  by  SAS  Institute  Inc.,  is  used  to  generate 
the  required  runs  and  perform  the  subsequent  analysis  for  the  2®r4  fractional  factorial 
experiment.  Using  the  parameter  ranges  defined  in  Sec.  3.3.1,  JMP  generates  the 
runs  shown  in  Tab.  3.4. 


Once  the  runs  are  defined  by  JMP,  Xpatch®  is  used  to  generate  ISAR  images 
(stored  as  binary  .trace  hies)  for  each  geometry  and  set  of  runs.  The  pixel  data  is 
extracted  from  the  ISAR  images  using  the  si2rgb.x  utility.  The  resulting  .iq  hie 
(the  pixel  magnitude  of  the  ISAR  image)  is  then  read  into  Matlab®  and  analyzed 
for  cross  range  smear  where  each  cross  range  slice  is  analyzed  for  cross  range  smear 
per  Sec.  2.5.  The  maximum  smear  for  each  image  is  the  smear  value  used  for  the 
subsequent  DOE  analysis  and  input  parameter  characterization. 


The  same  statistical  program,  JMP,  is  then  used  to  construct  a  mathematical 
model  of  the  SBR  and  ISAR  imaging  process.  Given  all  second-order  parameters 
are  not  included  in  the  model  due  to  DOF  limitations,  the  included  second-order 
parameters  are  varied  to  produce  a  model  of  least  squares  best  ht.  The  ht  of  each 
model  is  described  by  its  R2  value,  where  an  R 2  approaching  unity  indicates  a  perfect 
ht  with  empirical  data,  i.e.,  the  calculated  parameter  coefficients  accurately  charac¬ 
terize  cross  range  smear.  More  formally,  R2  is  the  proportion  of  the  experimental 
system  described  by  the  DOE-calculated  model,  defined  as  the  sum  of  squared  due 
to  regression  (SSR)  divided  by  the  total  sum  of  squares  (SST)  [9].  The  R2  measure 
is  written  as 


El=1  (smearmodei,n  ~  coY 


R2  = 


(3.15) 


J2u=l  (smear actual, n  -  c0)2 
where  N  is  the  number  of  runs  accomplished  and  Co  is  the  constant  parameter  in 
Eqn.  (3.4). 
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Table  3.4:  A  list  of  the  required  runs  for  the  2y  4  fractional  factorial  experiment  design,  generated  with  JMP.  The 
Pattern  column  identifies  high  (+),  low  (-),  and  center  (0)  parameter  values. 


Run 

DOE 

Pattern 

Ray 

Density 

(rays/A) 

Frequency 

Extent 

(GHz) 

Frequency 

Center 

(GHz) 

Frequency 

Increment 

(GHz) 

Angular 

Span 

(Degrees) 

Angular 

Increment 

(Degrees) 

First  Bounce 
Algorithm 

(PO/SBR) 

Diffraction 
Edge  File 

(Yes/No) 

1 

5 

1 

10 

0.005 

5 

0.01 

PO 

No 

2 

_ +  +  +  +- 

5 

1 

10 

0.01 

20 

0.1 

SBR 

No 

3 

- +  -  +  +  -  + 

5 

1 

15 

0.005 

20 

0.1 

PO 

Yes 

4 

- +  + - ++ 

5 

1 

15 

0.01 

5 

0.01 

SBR 

Yes 

5 

-  + - +  -  ++ 

5 

10 

10 

0.005 

20 

0.01 

SBR 

Yes 

6 

-  +  -  +  -  +  -  + 

5 

10 

10 

0.01 

5 

0.1 

PO 

Yes 

7 

-  +  + - +  +- 

5 

10 

15 

0.005 

5 

0.1 

SBR 

No 

8 

_++++___ 

5 

10 

15 

0.01 

20 

0.01 

PO 

No 

9 

+ _ +  ++ 

20 

1 

10 

0.005 

5 

0.1 

SBR 

Yes 

10 

+ - +  + - + 

20 

1 

10 

0.01 

20 

0.01 

PO 

Yes 

11 

+  -  +  -  +  -+- 

20 

1 

15 

0.005 

20 

0.01 

SBR 

No 

12 

+  -  +  +  -  + - 

20 

1 

15 

0.01 

5 

0.1 

PO 

No 

13 

+  + - +  + - 

20 

10 

10 

0.005 

20 

0.1 

PO 

No 

14 

+  +  -  + - +- 

20 

10 

10 

0.01 

5 

0.01 

SBR 

No 

15 

+  +  + _ + 

20 

10 

15 

0.005 

5 

0.01 

PO 

Yes 

16 

++++++++ 

20 

10 

15 

0.01 

20 

0.1 

SBR 

Yes 

17 

0  0  0  0  0  0  - 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

No 

18 

0  0  0  0  0  0  -  + 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

Yes 

19 

ooooooi — 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

No 

20 

0  0  0  0  0  0  ++ 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

Yes 

3.5  Cross  Range  Smear  Reduction 

In  addition  to  the  input  parameter  characterization  described  in  Sec.  3.4,  an 
examination  of  the  smear  reduction  method  proposed  by  Bhalla  and  Ling  [2]  is 
performed.  The  proposed  smear  reduction  implements  a  sliding  window  filter  across 
constant  frequency  arcs  of  the  spatial  frequency  annulus.  This  aspect  averaging 
reduces  discontinuities  in  the  spatial  frequency  data  that,  according  to  Bhalla  and 
Ling,  are  the  angular  scintillations  that  cause  cross  range  smear  in  the  resulting 
ISAR  image. 

The  averaging  technique  calls  for  a  3-aspect  sliding  window  applied  across 
each  spatial  frequency  aspect  arc.  Before  the  windowing  can  be  applied,  however, 
the  spatial  frequency  data  is  exported  to  a  text  (.field)  file  using  ss2field.x.  The 
data  is  then  read  into  Matlab®  and  windowed.  New  data  is  written  out  as  another 
.field  file,  where  it  is  reinserted  into  the  Xpatch®  SBR/ISAR  process  (illustrated  in 
Fig.  2.4)  using  the  ss2trace.x  utility  to  generate  the  output  ISAR  image.  This 
command  is  the  same  as  used  to  generate  the  ISAR  image  from  the  unfiltered  data. 
Cross  range  smear  is  then  analyzed  as  per  Sec.  3.4.  Cross  range  smear  in  these  new 
images  is  compared  to  those  that  have  not  been  processed  with  this  smear  reduction 
technique. 

In  addition  to  the  proposed  3-aspect  sliding  window,  a  9-aspect  sliding  window 
is  implemented  for  comparison.  This  extended  window  is  implemented  in  the  same 
manner  as  the  3-aspect  window,  and  the  analysis  is  performed  identically. 


54 


IV.  Cross  Range  Smear  Characterization  Results 

4-1  Introduction 

This  chapter  presents  the  results  and  analysis  of  cross  range  smear  character¬ 
ization  in  Xpatch®  Inverse  Synthetic  Aperture  Radar  (ISAR)  images.  These  ISAR 
images  are  generated  from  the  geometries  and  input  parameters  of  Chap.  III.  For 
each  geometry,  a  predesignated  set  of  input  parameter  values  is  used  to  generate 
ISAR  images.  These  images  are  then  individually  analyzed  for  cross  range  smear. 
These  smear  values  serve  as  the  basis  for  the  empirical  smear  model  for  that  geom¬ 
etry.  Due  to  the  high  dependance  of  geometry  on  electromagnetic  (EM)  scattering, 
each  geometry  has  its  own  empirical  smear  model  and  is  individually  discussed. 

For  each  geometry,  an  empirical  cross  range  smear  model  is  presented  and 
discussed.  This  model,  generated  with  Design  of  Experiments  (DOE)  techniques, 
approximates  the  smear  computed  in  an  Xpatch®  ISAR  image  given  a  set  of  param¬ 
eter  inputs.  Each  model’s  parameter  coefficients  characterize  how  those  parameters 
affect  cross  range  smear  given  that  scattering  geometry. 

A  separate  cross  range  smear  model  is  necessary  for  each  geometry  since  EM 
scattering,  and  therefore  ISAR  images,  are  highly  dependent  on  shape.  The  pa¬ 
rameters  included  in  each  model  are  varied  to  best  fit  the  empirical  model  to  the 
experimental  smear  data.  The  parameters  and  their  coefficients  form  the  smear 
model  as 

smear  =  f{pi,P2,P3,P4,P5,Pa,P7,Ps)  (4.1) 

=  Co  +  C1P1  +  C-2P-2  +  C3p3  +  c4p4  +  C5p5  +  c6p6  +  c7p7  +  C-sPs 
+Cupl  +  C12P1P2  +  CV3PlP3  +  C14P1PA 

+C22P2  +  C23P2P3  +  C24P2P4  +  •  •  •  (4-2) 
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where  the  actual  included  first  and  second-order  parameters  and  coefficients  vary 
with  each  geometries’  model.  It  is  important  to  note  that  parameters  are  linearly 
mapped  to  vary  between  —1  and  1,  so  parameter  coefficients  are  directly  comparable. 
The  constant  Co  is  referred  to  as  the  intercept,  and  is  included  in  each  DOE  model 
for  functional  completeness. 

All  first-order  parameters  are  included  in  each  model  to  serve  as  a  basis  for 
comparison,  but  additional  parameters  are  changed  to  enhance  R 2  fit  and  illustrate 
effect. 

4.2  Flat  Rectangular  Plate 

The  flat  rectangular  plate  produces  ISAR  images  with  smear  values  given  in 
Tab.  4.1.  Using  this  data,  an  empirical  model  with  the  parameters  and  parameter 
coefficients  shown  in  Tab.  4.2  is  generated.  The  model  is  very  well  fit  due  as  indicated 
by  its  high  R2  value. 

The  amount  of  smear  varies  between  between  1500  and  3000.  The  intercept  is 
by  far  the  largest  term,  and  parameter  coefficients  are  small  in  comparison.  How¬ 
ever,  parameter  coefficients  are  still  distinguishable  and  meaningfully  contribute  to 
the  smear  model.  The  largest  parameter  contribution,  by  far,  is  frequency  center. 
Bandwidth  is  just  over  20%  less  important,  followed  by  the  squared  ray  density  term. 

The  ray  density  squared  contribution  is  the  third  strongest,  meaning  the  model 
has  significant  curvature  over  its  ray  density  dimension.  The  second-order  ray  density 
parameter,  when  combined  with  angular  increment,  frequency  increment,  and  angu¬ 
lar  span,  are  also  of  note,  although  less  so  than  the  ray  density  squared  term.  The 
first-order  ray  density  term  is  almost  insignificant.  Besides  the  ray  density  squared 
term,  therefore,  ray  density  does  not  contribute  a  large  amount  to  the  smear  model. 
This  ray  density  unimportance  is  attributable  to  the  perfectly  flat  geometry,  where 
ray  discretization  effects  are  reduced. 
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Table  4.1:  The  smear  value  for  each  flat  rectangular  plate  ISAR  image  is  given  with  its  constituent  input  parameters. 
Each  ISAR  image  is  generated  using  Xpatch®  with  the  flat  rectangular  plate  geometry,  described  in  Sec.  3.2.1.  The 
Pattern  column  identifies  high  (+),  low  (-),  and  center  (0)  parameter  values. 


Run 

DOE 

Pattern 

Ray 

Density 

(rays/A) 

Frequency 

Extent 

(GHz) 

Frequency 

Center 

(GHz) 

Frequency 

Increment 

(GHz) 

Angular 

Span 

(Degrees) 

Angular 

Increment 

(Degrees) 

First  Bounce 
Algorithm 

(PO/SBR) 

Diffraction 
Edge  File 

(Yes/No) 

Smear 

1 

5 

1 

10 

0.005 

5 

0.01 

PO 

No 

1903.02 

2 

_ +  +  +  +- 

5 

1 

10 

0.01 

20 

0.1 

SBR 

No 

1619.73 

3 

- +  -  +  +  -  + 

5 

1 

15 

0.005 

20 

0.1 

PO 

Yes 

2671.94 

4 

- +  + - ++ 

5 

1 

15 

0.01 

5 

0.01 

SBR 

Yes 

2312.11 

5 

-  + - +  -  ++ 

5 

10 

10 

0.005 

20 

0.01 

SBR 

Yes 

2591.64 

6 

-  +  -  +  -  +  -  + 

5 

10 

10 

0.01 

5 

0.1 

PO 

Yes 

2145.23 

7 

-  +  + - +  +- 

5 

10 

15 

0.005 

5 

0.1 

SBR 

No 

3553.98 

8 

_++++___ 

5 

10 

15 

0.01 

20 

0.01 

PO 

No 

3830.17 

9 

+ _ +  ++ 

20 

1 

10 

0.005 

5 

0.1 

SBR 

Yes 

1736.79 

10 

+ - +  + - + 

20 

1 

10 

0.01 

20 

0.01 

PO 

Yes 

1551.13 

11 

+  -  +  -  +  -+- 

20 

1 

15 

0.005 

20 

0.01 

SBR 

No 

2960.91 

12 

+  -++-  +  -- 

20 

1 

15 

0.01 

5 

0.1 

PO 

No 

2672.02 

13 

+  + - +  + - 

20 

10 

10 

0.005 

20 

0.1 

PO 

No 

2671.73 

14 

+  +  -  + - +- 

20 

10 

10 

0.01 

5 

0.01 

SBR 

No 

2163.46 

15 

+  +  + _ + 

20 

10 

15 

0.005 

5 

0.01 

PO 

Yes 

3197.39 

16 

++++++++ 

20 

10 

15 

0.01 

20 

0.1 

SBR 

Yes 

4393.26 

17 

0  0  0  0  0  0  - 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

No 

2894.72 

18 

0  0  0  0  0  0  — b 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

Yes 

2913.16 

19 

0  0  0  0  0  0  +- 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

No 

2956.15 

20 

0  0  0  0  0  0  ++ 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

Yes 

2906.84 

Table  4.2:  Flat  Rectangular  Plate  Smear  Model 


Parameter 

Coefficient 

Intercept 

2917.72 

Frequency  Center 

575.57 

Bandwidth 

444.95 

Ray  Density*Ray  Density 

-294.31 

Angular  Span 

162.91 

Ray  Density*Angular  Increment 

140.44 

Ray  Density*Frequency  Increment 

64.15 

Ray  Density* Angular  Span 

63.01 

Frequency  Increment* Angular  Increment 

61.99 

Angular  Increment 

59.68 

Ray  Density 

44.93 

Diffraction  Edge  File 

-40.32 

Frequency  Increment 

-37.52 

First  Bounce  Algorithm 

37.22 

R2 

0.9670 

For  the  flat  rectangular  plate,  the  first  bounce  algorithm,  diffraction  edge  file, 
and  frequency  increment  are  all  relatively  insignificant.  However,  since  they  are  both 
first-order  parameters,  they  are  included  as  a  comparison  against  other  parameter 
effects. 

While  the  flat  rectangular  plate  may  not  be  the  most  interesting  target,  it 
serves  as  an  excellent  comparison  with  the  ridged  plate  of  Sec.  4.3. 


4-3  Ridged  Plate 

The  ridged  plate  produces  ISAR  images  with  smear  values  given  in  Tab.  4.3. 
Using  this  data,  an  empirical  model  with  the  parameters  and  parameter  coefficients 
shown  in  Tab.  4.4  is  generated.  The  model  is  very  well  fit,  evidenced  by  its  high  R 2 
value. 

The  amount  of  smear  varies  significantly  between  runs;  the  standard  deviation 
over  all  the  smear  values  is  just  over  1400.  This  variation  indicates  the  presence  of 
strong  parameter  coefficients,  agreeing  with  the  large  numbers  seen  in  Tab.  4.4. 
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Table  4.3:  The  smear  value  for  each  ridged  plate  ISAR  image  is  given  with  its  constituent  input  parameters.  Each 
ISAR  image  is  generated  using  Xpatch®  with  the  ridged  plate  geometry,  described  in  Sec.  3.2.2.  The  Pattern  column 
identifies  high  (+),  low  (-),  and  center  (0)  parameter  values. 


Run 

DOE 

Pattern 

Ray 

Density 

(rays/A) 

Frequency 

Extent 

(GHz) 

Frequency 

Center 

(GHz) 

Frequency 

Increment 

(GHz) 

Angular 

Span 

(Degrees) 

Angular 

Increment 

(Degrees) 

First  Bounce 
Algorithm 

(PO/SBR) 

Diffraction 
Edge  File 

(Yes/No) 

Smear 

1 

5 

1 

10 

0.005 

5 

0.01 

PO 

No 

435.56 

2 

_ +  +  +  +- 

5 

1 

10 

0.01 

20 

0.1 

SBR 

No 

419.96 

3 

- +  -  +  +  -  + 

5 

1 

15 

0.005 

20 

0.1 

PO 

Yes 

294.64 

4 

- +  + - ++ 

5 

1 

15 

0.01 

5 

0.01 

SBR 

Yes 

293.59 

5 

-  + - +  -  ++ 

5 

10 

10 

0.005 

20 

0.01 

SBR 

Yes 

2885.33 

6 

-  +  -  +  -  +  -  + 

5 

10 

10 

0.01 

5 

0.1 

PO 

Yes 

408.07 

7 

-  +  + - +  +- 

5 

10 

15 

0.005 

5 

0.1 

SBR 

No 

641.98 

8 

_++++___ 

5 

10 

15 

0.01 

20 

0.01 

PO 

No 

4501.35 

9 

+ _ +  ++ 

20 

1 

10 

0.005 

5 

0.1 

SBR 

Yes 

380.54 

10 

+ - +  + - + 

20 

1 

10 

0.01 

20 

0.01 

PO 

Yes 

393.46 

11 

+  -  +  -  +  -+- 

20 

1 

15 

0.005 

20 

0.01 

SBR 

No 

274.16 

12 

+  -++-  +  -- 

20 

1 

15 

0.01 

5 

0.1 

PO 

No 

248.76 

13 

+  + - +  + - 

20 

10 

10 

0.005 

20 

0.1 

PO 

No 

3108.54 

14 

+  +  -  + - +- 

20 

10 

10 

0.01 

5 

0.01 

SBR 

No 

378.75 

15 

+  +  + _ + 

20 

10 

15 

0.005 

5 

0.01 

PO 

Yes 

630.39 

16 

++++++++ 

20 

10 

15 

0.01 

20 

0.1 

SBR 

Yes 

4748.81 

17 

0  0  0  0  0  0  - 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

No 

1046.86 

18 

0  0  0  0  0  0  — b 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

Yes 

1066.90 

19 

0  0  0  0  0  0  +- 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

No 

1059.84 

20 

0  0  0  0  0  0  ++ 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

Yes 

1080.58 

Table  4.4:  Ridged  Plate  Smear  Model 


Parameter 

Coefficient 

Intercept 

1063.55 

Bandwidth 

910.16 

Angular  Span 

825.54 

Ray  Density*Angular  Increment 

822.57 

Ray  Density *Diffraction  Edge  File 

266.26 

Frequency  Center 

201.47 

Ray  Density*Ray  Density 

189.20 

Ray  Density*First  Bounce  Algorithm 

174.99 

Frequency  Increment 

171.35 

Angular  Increment 

28.67 

Ray  Density 

17.68 

Diffraction  Edge  File 

3.33 

First  Bounce  Algorithm 

1.45 

R2 

0.9989 

Similar  to  the  flat  rectangular  plate  of  Sec.  4.2,  the  dominant  model  term  is  the 
intercept.  The  strongest  parameter  contribution  is  bandwidth,  followed  by  angular 
span  and  then  ray  density* angular  increment.  Other  parameters  are  also  significant, 
but  contribute  considerably  less  than  these  first  three  parameters. 

Of  note,  ray  density  does  not  significantly  factor  into  the  model  as  a  first-order 
effect  but  does  have  a  very  large  contribution  when  combined  as  a  second-order 
parameter  with  diffraction  edge  file,  first  bounce  algorithm,  and  especially  angular 
increment.  The  ray  density  squared  term  is  also  significant,  contributing  curvature 
to  the  model  across  that  dimension.  Ray  density  therefore  does  not  have  a  large 
linear  contribution  to  the  smear  for  this  geometry,  but  becomes  an  important  factor 
when  considered  together  with  angular  increment. 

As  a  group,  the  frequency  dependant  terms  (bandwidth,  frequency  center,  and 
frequency  increment)  are  also  important  to  presence  of  smear.  Bandwidth,  in  this 
case,  is  by  far  the  largest  contributor  of  this  group.  Compared  to  the  flat  rectangular 
plate  parameter  coefficients,  frequency  center  is  much  less  important  in  this  model; 
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frequency  center  has  the  largest  parameter  coefficient  in  the  flat  rectangular  plate 
smear  model  but  is  relegated  to  fifth  largest  here. 

The  use  of  a  diffraction  edge  file  or  change  in  the  first  bounce  algorithm  has 
very  little  effect  on  the  presence  of  smear. 

4-4  Slicy 

The  slicy  geometry  produces  ISAR  images  with  smear  values  given  in  Tab.  4.5. 
Using  this  data,  an  empirical  model  with  the  parameters  and  parameter  coefficients 
shown  in  Tab.  4.6  is  generated.  The  model  is  very  well  fit,  evidenced  by  the  high  R 2 
value. 

The  amount  of  smear  generated  using  this  geometry  does  not  vary  much  be¬ 
tween  runs;  the  standard  deviation  over  all  smear  values  is  just  below  230.  This 
lack  of  variation  results  in  low  parameter  coefficients.  Table  4.6  agrees,  showing  that 
parameter  coefficients  are  all  much  smaller  than  the  intercept  term. 

By  far,  the  dominant  model  term  is  the  intercept.  The  strongest  parameter 
contribution  is  the  frequency  center,  followed  by  moderately  closely  by  the  ray  den¬ 
sity  squared  term  and  bandwidth.  The  strength  of  subsequent  parameters  decreases 
greatly  from  there. 

As  mentioned,  the  ray  density  squared  term  is  very  strong,  contributing  a 
lot  of  curvature  to  the  model’s  ray  density  dimension.  However,  the  lack  of  ray 
density’s  first-order  importance  means  that  over  the  measured  span  (5  to  20  rays  per 
wavelength)  the  difference  between  the  model’s  low  and  high  side  ray  density  response 
is  relatively  small.  However,  the  ray  density  has  a  larger  effect  when  combined  with 
angular  increment  and  frequency  center  as  a  second-order  term  . 

Bandwidth  and  frequency  center  both  have  strong  contributions  to  the  extent 
of  cross  range  smear.  This  contribution  is  similar  to  other  models.  Below  these  two 
parameters  falls  angular  span,  also  consistent  across  models. 
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Table  4.5:  The  smear  value  for  each  slicy  ISAR  image  is  given  with  its  constituent  input  parameters.  Each  ISAR  image 
is  generated  using  Xpatch®  with  the  slicy  geometry,  described  in  Sec.  3.2.3.  The  Pattern  column  identifies  high  (+),  low 
(-),  and  center  (0)  parameter  values. 


Run 

DOE 

Pattern 

Ray 

Density 

(rays/A) 

Frequency 

Extent 

(GHz) 

Frequency 

Center 

(GHz) 

Frequency 

Increment 

(GHz) 

Angular 

Span 

(Degrees) 

Angular 

Increment 

(Degrees) 

First  Bounce 
Algorithm 

(PO/SBR) 

Diffraction 
Edge  File 

(Yes/No) 

Smear 

1 

5 

1 

10 

0.005 

5 

0.01 

PO 

No 

481.08 

2 

_ +  +  +  +- 

5 

1 

10 

0.01 

20 

0.1 

SBR 

No 

489.66 

3 

- +  -  +  +  -  + 

5 

1 

15 

0.005 

20 

0.1 

PO 

Yes 

821.29 

4 

- +  + - ++ 

5 

1 

15 

0.01 

5 

0.01 

SBR 

Yes 

743.10 

5 

-  + - +  -  ++ 

5 

10 

10 

0.005 

20 

0.01 

SBR 

Yes 

766.74 

6 

-  +  -  +  -  +  -  + 

5 

10 

10 

0.01 

5 

0.1 

PO 

Yes 

553.24 

7 

-  +  + - +  +- 

5 

10 

15 

0.005 

5 

0.1 

SBR 

No 

1041.39 

8 

_++++___ 

5 

10 

15 

0.01 

20 

0.01 

PO 

No 

1113.69 

9 

+ _ +  ++ 

20 

1 

10 

0.005 

5 

0.1 

SBR 

Yes 

512.72 

10 

+ - +  + - + 

20 

1 

10 

0.01 

20 

0.01 

PO 

Yes 

470.47 

11 

+  -  +  -  +  -+- 

20 

1 

15 

0.005 

20 

0.01 

SBR 

No 

696.93 

12 

+  -++-  +  -- 

20 

1 

15 

0.01 

5 

0.1 

PO 

No 

812.63 

13 

+  + - +  + - 

20 

10 

10 

0.005 

20 

0.1 

PO 

No 

822.94 

14 

+  +  -  + - +- 

20 

10 

10 

0.01 

5 

0.01 

SBR 

No 

534.72 

15 

+  +  + _ + 

20 

10 

15 

0.005 

5 

0.01 

PO 

Yes 

990.57 

16 

++++++++ 

20 

10 

15 

0.01 

20 

0.1 

SBR 

Yes 

1286.60 

17 

0  0  0  0  0  0  - 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

No 

914.87 

18 

0  0  0  0  0  0  — b 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

PO 

Yes 

873.20 

19 

0  0  0  0  0  0  +- 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

No 

899.67 

20 

0  0  0  0  0  0  ++ 

12.5 

5.5 

12.5 

0.0075 

12.5 

0.055 

SBR 

Yes 

883.85 

Table  4.6:  Slicy  Smear  Model 


Parameter 

Coefficient 

Intercept 

892.90 

Frequency  Center 

179.66 

Ray  Density*Ray  Density 

-134.29 

Bandwidth 

130.13 

Ray  Density* Angular  Increment 

52.36 

Angular  Span 

49.93 

Bandwidth*Frequency  Center 

39.66 

Angular  Increment 

33.95 

Ray  Density*Frequency  Increment 

18.25 

Angular  Span*  Angular  Increment 

12.63 

Frequency  Increment 

-8.10 

Ray  Density 

7.34 

First  Bounce  Algorithm*Diffraction  Edge  File 

6.46 

Diffraction  Edge  File 

4.71 

Ray  Density* Frequency  Center 

1.07 

First  Bounce  Algorithm 

0.07 

R2 

0.9968 

The  first-order  effects  of  a  diffraction  edge  file  or  first  bounce  algorithm  change 
is  next  to  negligible.  This  (lack  of)  effect  is  consistent  with  the  results  from  other 
geometries’  cross  range  smear  models. 


4-5  Conclusions 

Certain  factors,  such  as  target  geometry  and  a  few  input  parameters,  influence 
cross  range  smear  a  great  deal.  Most  factors,  however,  have  a  varying  effect  on 
smear.  The  choice  of  Xpatch®’s  first  bounce  algorithm  or  diffraction  edge  hie  has 
almost  no  effect  on  cross  range  smear. 

The  primary  conclusion  regarding  cross  range  smear  is  that,  as  expected,  cross 
range  smear  is  primarily  a  function  of  geometry.  This  conclusion  is  supported  by  the 
large  difference  in  smear  values  between  geometries.  In  addition,  since  ISAR  images 


63 


vary  depending  on  the  scattering  geometry1,  comparing  the  amount  of  cross  range 
smear  in  images  from  one  geometry  to  a  dramatically  different  one  is  of  limited  value. 
However,  the  variation  of  input  parameters  while  keeping  the  geometry  constant 
creates,  through  DOE  analysis,  smear  characteristics  comparable  across  geometries. 

Ray  density,  as  a  first-order  effect,  is  consistently  a  minor  influence  on  cross 
range  smear  regardless  of  geometry.  However,  when  considered  as  a  second-order 
term  with  parameters  such  as  angular  increment,  it  has  a  much  larger  effect.  In  ad¬ 
dition,  there  is  consistently  model  curvature  induced  by  a  strong  ray  density  squared 
term.  These  effects  become  more  influential  as  the  target  complexity  increases  be¬ 
yond  the  flat  rectangular  plate.  This  result  is  consistent  with  Bhalla  and  Ling’s  [2] 
explanation  of  cross  range  smear  cause  as  angular  scintillation. 

The  bandwidth  and  frequency  center  terms  are  also  important  to  the  creation 
of  cross  range  smear.  In  Xpatch®,  ray  tubes  are  bounced  from  each  aspect  and 
the  scattered  held  as  a  function  of  frequency  is  known  in  closed-form  once  the  ray 
bounces  are  complete.  There  should  therefore  not  be  any  loss  of  coherency  or  un¬ 
expected  discontinuity  over  frequency  at  a  constant  aspect.  However,  a  general 
increase  in  frequency  (which  both  bandwidth  and  frequency  center  determine)  will 
increase  the  scattering  target’s  electrical  size,  meaning  that  an  effective  scattered 
held  representation  requires  a  finer  target  surface  sampling.  This  requirement  may 
induce  unexpected  but  otherwise  minimized  aspect-based  discontinuities,  and  is  not 
specifically  addressed  in  Bhalla  and  Ling’s  [2]  explanation. 

Cross  range  smear  is  characterized  in  terms  of  input  parameters  for  a  limited 
number  of  geometries,  so  expanding  this  characterization  to  cover  several  geometries 
needs  to  be  done  carefully.  Besides  the  natural  dependance  on  geometry,  frequency 
center  and  bandwidth  exert  an  important  influence  on  cross  range  smear  presence. 

lrThe  variation  of  ISAR  images  depending  on  target  geometry  makes  it  a  powerful  scattering 
signature. 
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Ray  density  is  also  a  very  important  factor,  but  only  as  a  combined,  second-order 
parameter. 
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V.  Cross  Range  Smear  Reduction  Results 

5. 1  Introduction 

This  chapter  presents  the  results  and  analysis  of  the  smear  reduction  technique 
proposed  by  Bhalla  and  Ling  [2]  and  described  in  Sec.  3.5.  A  smear  reduction 
analysis  is  presented  for  each  geometry  used  in  smear  characterization,  and  a  general 
discussion  applicable  to  all  geometries  follows.  In  addition  to  the  proposed  smear 
reduction  technique,  a  similar  technique  using  a  9-aspect  extended  sliding  window 
is  also  analyzed. 

The  smear  reduction  technique  is  applied  to  each  Xpatch®-generated  scattered 
field  file,  and  the  result  used  by  Xpatch®  to  produce  an  ISAR  image.  Each  image 
is  evaluated  for  cross  range  smear  as  per  Sec.  2.5,  and  the  results  compared  to  the 
unfiltered  ISAR  image’s  cross  range  smear  values.  For  consistency,  each  ISAR  image 
is  shown  with  a  50  dB  dynamic  range. 

5.2  Flat  Rectangular  Plate 

The  flat  rectangular  plate  displays  a  fair  amount  of  cross  range  smear,  and 
the  smear  reduction  technique  (as  originally  proposed)  performs  well.  Both  the 
maximum  cross  range  smear,  assigned  as  the  smear  for  the  entire  image,  and  the 
mean  cross  range  smear  are  reduced  in  all  ISAR  images/runs.  Figure  5.1  shows  the 
maximum  cross  range  smear  over  each  run,  while  Fig.  5.2  shows  the  corresponding 
mean  values. 

Using  the  original  3-aspect  smear  reduction  technique,  maximum  smear  is  re¬ 
duced  in  every  ISAR  image  except  two.  The  amount  of  smear  reduction  varies, 
however,  and  is  not  a  linear  function  of  smear  value.  The  9-aspect  smear  reduction 
reduced  smear  even  further. 
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Figure  5.1:  A  comparison  of  flat  rectangular  plate  maximum 
cross  range  smear  reduction  by  run  number.  The  maximum  cross 
range  smear  is  the  smear  value  assigned  to  the  entire  image  for  the 
smear  characterization. 


run 

Figure  5.2:  A  comparison  of  flat  rectangular  plate  mean  cross 
range  smear  reduction  by  run  number.  The  mean  cross  range  smear 
is  computed  over  each  ISAR  image’s  entire  range. 
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original  and  3-aspect  filtered  ISAR  im¬ 
ages  for  flat  rectangular  plate  run  16. 
Given  this  geometry,  cross  range  smear 
is  reduced  by  the  largest  amount  in  this 
run. 


Figure  5.4:  The  difference  in  cross 
range  smear  between  the  3-aspect  fil¬ 
tered  and  original  flat  rectangular  plate 
run  16  ISAR  images.  This  plot  is  the  dif¬ 
ference  between  the  two  smear  functions 
in  Fig.  5.3. 


The  smear  reduction  technique  performs  best  in  run  16,  which  also  has  the 
largest  (maximum)  smear  value.  The  original  and  3-aspect  filtered  ISAR  images 
are  shown  in  Fig.  5.7  and  5.8,  respectively.  Detailed  in  Fig.  5.3,  smear  is  almost 
negligible  through  all  range  cells  except  for  where  the  flat  plate  scattering  return  is 
located.  Small  smear  sidelobes  are  present  due  to  smeared  range-based  sidclobes  in 
the  ISAR  images,  but  are  much  weaker  than  the  peak  smear  value.  However,  smear 
reduction  is  fairly  constant  across  the  entire  range,  through  the  peak  and  over  the 
sidelobes.  This  consistency  is  illustrated  by  Fig.  5.4,  and  visually  noticeable  in  the 
filtered  ISAR  image’s  reduced  sidelobes. 

The  9-aspect  filtered  ISAR  image  is  shown  in  Fig.  5.9.  Notably,  the  smear 
present  above  and  below  the  unfiltered  plate  is  reduced  below  the  image’s  dynamic 
range.  Compared  to  the  original,  maximum  smear  is  reduced  81%  and  mean  smear 
reduced  by  80%.  The  difference  in  smear  is  illustrated  in  Figure  5.5  and  5.6. 

The  smallest  smear  reduction,  or  rather  the  largest  increase,  occurs  in  the  image 
generated  with  run  5.  The  original  ISAR  image  is  shown  in  Fig.  5.14.  After  filtering 
with  a  3-aspect  sliding  window,  the  ISAR  image  shown  in  Fig.  5.15  is  produced. 
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Figure  5.5:  Cross  range  smear  in  the 
original  and  9-aspect  filtered  ISAR  im¬ 
ages  for  flat  rectangular  plate  run  16. 
Given  this  geometry,  cross  range  smear 
is  reduced  by  the  largest  amount  in  this 
run. 


Figure  5.6:  The  difference  in  cross 
range  smear  between  the  9-aspect  fil¬ 
tered  and  original  flat  rectangular  plate 
run  16  ISAR  images.  This  plot  is  the  dif¬ 
ference  between  the  two  smear  functions 
in  Fig.  5.5. 
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Figure  5.7:  The  flat  rectangular  plate  ISAR  image  generated 
with  the  parameters  given  for  run  16.  This  image  is  unfiltered. 
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Figure  5.8:  The  flat  rectangular  plate  ISAR  image  generated 
with  the  parameters  given  for  run  16.  This  image  is  filtered  with 
a  3-aspect  sliding  window,  as  originally  proposed.  Maximum 
cross  range  smear  is  reduced  by  47%. 
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Figure  5.9:  The  flat  rectangular  plate  ISAR  image  generated 
with  the  parameters  given  for  run  16.  This  image  is  filtered 
with  a  9-aspect  sliding  window,  reducing  maximum  cross  range 
smear  by  81%. 
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Figure  5.10:  Cross  range  smear  in  the 
original  and  3-aspect  filtered  ISAR  im¬ 
ages  for  flat  rectangular  plate  run  5. 
Given  this  geometry,  cross  range  smear 
is  reduced  by  the  smallest  amount  in  this 
run.  Maximum  smear  actually  increased 
by  a  small  amount,  but  mean  smear  de¬ 
creased  by  a  similarly  small  margin. 
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Figure  5.11:  The  difference  in  cross 
range  smear  between  the  3-aspect  fil¬ 
tered  and  original  flat  rectangular  plate 
run  5  ISAR  images.  This  plot  is  the  dif¬ 
ference  between  the  two  smear  functions 
in  Fig.  5.10. 


Given  the  maximum  smear  magnitude,  approximately  3800,  this  increase  is  all  but 
negligible  at  0.4%.  This  small  change  is  reinforced  by  Fig.  5.10,  which  shows  the 
amount  of  change  is  very  small  compared  to  the  magnitude  of  smear.  This  lack  of 
significant  change  is  seen  by  the  largely  overlapping  plot.  Figure  5.11  shows  the  net 
difference  between  the  two  smear  functions. 

When  filtered  with  a  9-aspect  sliding  window,  the  ISAR  image  in  Fig.  5.16  is 
produced.  Maximum  smear  is  reduced  from  the  original  image  by  2%,  contrasted 
with  the  slight  increase  using  the  3-aspect  window.  The  smear  difference  between 
this  image  and  the  original  is  detailed  in  Fig.  5.12  and  5.13.  This  lack  of  significant 
smear  reduction  results  in  nearly  identical  ISAR  images. 


5.3  Ridged  Plate 

The  amount  of  maximum  and  mean  smear  present  in  the  ridged  plate  runs 
varies  significantly,  as  did  the  smear  reduction  technique’s  performance.  In  general, 
smear  reduction  performs  best  on  runs  with  large  amounts  of  smear.  Also  of  note,  the 
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Figure  5.12:  Cross  range  smear  in  the 
original  and  9-aspect  filtered  ISAR  im¬ 
ages  for  flat  rectangular  plate  run  5. 
Given  this  geometry,  cross  range  smear 
is  reduced  by  the  smallest  amount  in  this 
run.  Maximum  smear  actually  increased 
by  a  small  amount,  but  mean  smear  de¬ 
creased  by  a  similarly  small  margin. 
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Figure  5.13:  The  difference  in  cross 
range  smear  between  the  9-aspect  fil¬ 
tered  and  original  flat  rectangular  plate 
run  5  ISAR  images.  This  plot  is  the  dif¬ 
ference  between  the  two  smear  functions 
in  Fig.  5.12. 
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Figure  5.14:  The  flat  rectangular  plate  ISAR  image  gener¬ 
ated  with  the  parameters  given  for  run  5.  This  image  is  unfil¬ 
tered. 
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Figure  5.15:  The  flat  rectangular  plate  ISAR  image  gener¬ 
ated  with  the  parameters  given  for  run  5.  This  image  is  filtered 
with  a  3-aspect  sliding  window,  as  originally  proposed.  Maxi¬ 
mum  cross  range  smear  is  increased  by  approximately  0.4%. 
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Figure  5.16:  The  flat  rectangular  plate  ISAR  image  gen¬ 
erated  with  the  parameters  given  for  run  5.  This  image  is 
filtered  with  a  9-aspect  sliding  window,  decreasing  maximum 
cross  range  smear  by  3%.  However,  with  this  window  size,  a 
new  low  level  constant  smear  is  induced.  This  new  smear  indi¬ 
cates  that  the  window  is  too  large  for  this  image. 
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Figure  5.17:  A  comparison  of  ridged  plate  maximum  cross  range 
smear  reduction  by  run  number.  The  maximum  cross  range  smear 
is  the  smear  value  assigned  to  the  entire  image  for  the  smear  char¬ 
acterization. 

mean  smear  values  are  reduced  more  consistently  than  the  maximum  smear  values. 
This  means  that  the  smear  reduction  technique  is  more  effective  at  reducing  off-peak 
levels  of  smear.  The  amount  of  smear  reduction,  in  either  the  maximum  or  mean 
case,  is  not  a  linear  function  of  smear  value.  Figure  5.17  shows  the  maximum  cross 
range  smear  over  each  run,  while  Fig.  5.18  shows  the  corresponding  mean  values. 

A  good  example  of  smear  reduction  is  run  19,  where  the  original  image  is  shown 
in  Fig.  5.23  and  the  3-aspect  filtered  image  in  Fig.  5.24.  The  two  ISAR  images 
produce  smear  functions  shown  in  Fig.  5.19.  This  run  does  not  have  the  largest 
maximum  or  mean  smear  reduction,  but  the  relatively  low  maximum  value  and 
moderate  mean  value  readily  show  smear  reduction  performance.  The  peak  value  is 
present  at  the  ridged  plate  location,  at  the  downrange  center,  but  a  dispersed  return 
further  downrange  is  also  smeared.  This  dispersed  response  is  due  to  the  ridged 
plate  geometry.  The  smear  reduction  technique  performs  well  for  the  smear  present 
on  both  the  maximum  and  dispersed  returns;  maximum  smear  is  reduced  by  9%. 
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Figure  5.18:  A  comparison  of  ridged  plate  mean  cross  range 
smear  reduction  by  run  number.  The  mean  cross  range  smear  is 
computed  over  each  ISAR  image’s  entire  range. 

The  largest  reduction,  though,  is  over  the  dispersed  region  just  downrange  of  the 
ridged  plate.  There,  the  dispersed  scattering  response  is  significantly  compressed 
towards  the  cross  range  center.  Even  though  this  dispersed  region  is  not  the  area  of 
maximum  smear,  as  seen  in  Fig.  5.20  it  is  the  area  of  largest  smear  reduction. 

The  9-aspect  filtered  image,  shown  in  Fig.  5.25  has  even  less  smear.  A  compar¬ 
ison  between  the  original  and  this  filtered  ISAR  image’s  smear  is  given  in  Fig.  5.21 
and  5.22.  Here,  maximum  smear  is  reduced  by  just  under  55%  and  mean  smear 
by  67%.  The  dispersed  region  is  noticeably  reduced  and  compressed  towards  the 
middle,  even  more  than  the  3-aspect  filtered  image. 

The  smallest  smear  reduction  occurs  in  the  image  generated  with  run  8,  from 
the  original  and  3-aspect  filtered  ISAR  images  shown  in  Fig.  5.30  and  5.31.  Maximum 
and  mean  smear  are  both  reduced,  as  shown  in  Fig.  5.26  and  5.27,  but  not  to  the 
extent  seen  in  other  runs.  The  reduction  through  the  dispersed  downrange  response 
is  not  very  evident  either,  resulting  in  two  nearly  identical  ISAR  images 
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Figure  5.19:  Cross  range  smear  in  the 
original  and  3-aspect  filtered  ISAR  im¬ 
ages  for  ridged  plate  run  19.  Cross  range 
maximum  and  mean  smear  are  both  suc¬ 
cessfully  reduced  between  the  original 
and  3- aspect  filtered  images.  This  run 
does  not  have  the  largest  amount  of 
smear  reduction,  but  overall  smear  re¬ 
duction  is  more  readily  evident  than  for 
other  runs  due  to  its  relatively  low  peak 
value. 


Figure  5.21:  Cross  range  smear  in  the 
original  and  9-aspect  filtered  ISAR  im¬ 
ages  for  ridged  plate  run  19.  Cross  range 
maximum  and  mean  smear  are  both  suc¬ 
cessfully  reduced.  This  run  does  not  have 
the  largest  amount  of  smear  reduction, 
but  overall  smear  reduction  is  more  read¬ 
ily  evident  than  for  other  runs  due  to  its 
relatively  low  peak  value. 


Figure  5.20:  The  difference  in  cross 
range  smear  between  the  3-aspect  fil¬ 
tered  and  original  ridged  plate  run  19 
ISAR  images.  This  plot  is  the  differ¬ 
ence  between  the  two  smear  functions  in 
Fig.  5.19. 


Figure  5.22:  The  difference  in  cross 
range  smear  between  the  9-aspect  fil¬ 
tered  and  original  ridged  plate  run  19 
ISAR  images.  This  plot  is  the  differ¬ 
ence  between  the  two  smear  functions  in 
Fig.  5.21. 
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Figure  5.23:  The  ridged  plate  ISAR  image  generated  with 
the  parameters  given  for  run  19.  This  image  is  unfiltered. 
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Figure  5.24:  The  Ridged  plate  ISAR  image  generated  with 
the  parameters  given  for  run  19.  This  image  is  filtered  with  a  3- 
aspect  sliding  window,  decreasing  maximum  cross  range  smear 
by  approximately  10%. 
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Figure  5.25:  The  Ridged  plate  ISAR  image  generated  with 
the  parameters  given  for  run  19.  This  image  is  filtered  with  a  9- 
aspect  sliding  window,  decreasing  maximum  cross  range  smear 
by  approximately  55%. 

The  9-aspect  filtered  image,  shown  in  Fig.  5.32  has  less  maximum  and  mean 
smear  than  its  3-aspect  filtered  counterpart.  A  comparison  between  the  original  and 
this  filtered  ISAR  image’s  smear  is  given  in  Fig.  5.28  and  5.29.  Here,  maximum 
smear  is  reduced  by  3%  and  mean  smear  by  30%  from  the  original.  The  dispersed 
region  is  noticeably  reduced,  much  more  than  the  original  or  the  the  3-aspect  filtered 
version.  This  reduction  is  reflected  in  the  mean  smear  reduction  and  is  noticeable 
in  the  ISAR  image  as  a  general  reduction  and  compression  towards  the  center.  The 
9-aspect  filter  performs  well  for  this  geometry. 

5-4  Slicy 

The  amount  of  maximum  and  mean  smear  present  in  the  slicy  runs  varies  a 
fair  amount.  The  standard  deviation  of  maximum  smear  over  all  runs  is  just  over 
200  with  an  mean  value  of  approximately  800.  The  cross  range  smear  reduction 
technique’s  performance  as  proposed  varied  by  quite  a  bit,  not  a  linear  function  of 
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Figure  5.26:  Cross  range  smear  in  the 
original  and  3-aspect  filtered  ISAR  im¬ 
ages  for  ridged  plate  run  8.  Given 
this  geometry,  cross  range  smear  is  re¬ 
duced  by  the  smallest  amount  in  this 
run.  Smear  reduction  is  still  evident,  es¬ 
pecially  through  the  range  cells  just  after 
the  peak. 


original  and  9-aspect  filtered  ISAR  im¬ 
ages  for  ridged  plate  run  8.  Given 
this  geometry,  cross  range  smear  is  re¬ 
duced  by  the  smallest  amount  in  this 
run.  Smear  reduction  is  still  evident,  es¬ 
pecially  through  the  range  cells  just  after 
the  peak. 


Figure  5.27:  The  difference  in  cross 
range  smear  between  the  3-aspect  fil¬ 
tered  and  original  ridged  plate  run  15 
ISAR  images.  This  plot  is  the  differ¬ 
ence  between  the  two  smear  functions  in 
Fig.  5.26. 


Figure  5.29:  The  difference  in  cross 
range  smear  between  the  9-aspect  fil¬ 
tered  and  original  ridged  plate  run  8 
ISAR  images.  This  plot  is  the  differ¬ 
ence  between  the  two  smear  functions  in 
Fig.  5.28. 
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Figure  5.30:  The  ridged  plate  ISAR  image  generated  with 
the  parameters  given  for  run  8.  This  image  is  unfiltered. 
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Figure  5.31:  The  ridged  plate  ISAR  image  generated  with 
the  parameters  given  for  run  8.  This  image  is  filtered  with  a  3- 
aspect  sliding  window,  decreasing  maximum  cross  range  smear 
by  approximately  2%. 
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Figure  5.32:  The  ridged  plate  ISAR  image  generated  with 
the  parameters  given  for  run  8.  This  image  is  filtered  with  a  Si- 
aspect  sliding  window,  decreasing  maximum  cross  range  smear 
by  approximately  3%. 

smear  value.  In  all  but  one  case,  cross  range  smear  is  reduced.  Figure  5.33  shows  the 
maximum  cross  range  smear  over  each  run,  while  Fig.  5.34  shows  the  corresponding 
mean  values. 

A  good  example  of  smear  reduction  is  run  12,  which  produces  the  original  ISAR 
image  shown  in  Fig.  5.39  and  the  3-aspect  filtered  image  in  Fig.  5.40.  These  images 
produce  smear  functions  shown  in  Fig.  5.35  and  directly  compared  in  Fig.  5.36.  The 
reduction  in  maximum  smear  is  the  second  largest  at  45%,  but  run  12  has  the  largest 
reduction  in  mean  smear.  Smear  reduction  is  noticeable  on  the  cross  range  smear 
located  at  the  slicy  model’s  front  and  middle,  as  well  as  over  the  range-based  lobing 
present  throughout  the  entire  image. 

The  9-aspect  filtered  image,  shown  in  Fig.  5.41  has  noticeably  less  maximum 
smear  than  its  3-aspect  filtered  counterpart.  A  comparison  between  the  original  and 
this  filtered  ISAR  image’s  smear  is  given  in  Fig.  5.37  and  5.38.  Here,  maximum  and 
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Figure  5.33:  A  comparison  of  slicy  maximum  cross  range  smear 
reduction  by  run  number.  The  maximum  cross  range  smear  is  the 
smear  value  assigned  to  the  entire  image  for  the  smear  characteri¬ 
zation. 
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Figure  5.34:  A  comparison  of  slicy  mean  cross  range  smear  re¬ 
duction  by  run  number.  The  mean  cross  range  smear  is  computed 
over  each  ISAR  image’s  entire  range. 
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Figure  5.35:  Cross  range  smear  in  the 
original  and  3-aspect  filtered  ISAR  im¬ 
ages  for  slicy  run  12.  Cross  range  maxi¬ 
mum  and  mean  smear  are  both  success¬ 
fully  reduced.  This  run  does  not  have 
the  largest  amount  of  smear  reduction, 
but  overall  smear  reduction  is  more  read¬ 
ily  evident  than  for  other  runs  due  to  its 
relatively  low  peak  value. 


Figure  5.36:  The  difference  in  cross 
range  smear  between  the  3-aspect  fil¬ 
tered  and  original  slicy  run  12  ISAR  im¬ 
ages.  This  plot  is  the  difference  between 
the  two  smear  functions  in  Fig.  5.35. 


smear  are  reduced  by  approximately  80%.  The  9-aspect  filter  performs  better  than 
the  3-aspect  filter  in  this  case. 

The  smallest  smear  reduction  using  the  3-aspect  sliding  window,  actually  a 
small  increase  in  maximum  and  mean  smear,  occurs  in  the  image  generated  with  run 
1.  The  smear  in  the  original  and  filtered  ISAR  images,  Fig.  5.46  and  5.47,  is  shown 
in  Fig.  5.42.  The  difference  in  smear,  further  detailed  in  Fig.  5.43,  is  very  small. 
This  small  difference  makes  the  two  ISAR  images  appear  virtually  indistinguishable. 

The  9-aspect  filtered  image,  shown  in  Fig.  5.48  has  a  little  more  smear  than 
its  3-aspect  filtered  counterpart.  A  comparison  between  the  original  and  this  filtered 
ISAR  image’s  smear  is  given  in  Fig.  5.44  and  5.45.  Here,  maximum  smear  is  increased 
by  2%  but  mean  smear  is  reduced  by  1%.  The  9-aspect  filter  performs  about  the 
same  as  the  3-aspect  filter  in  this  case. 
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Figure  5.37:  Cross  range  smear  in  the 
original  and  9-aspect  filtered  ISAR  im¬ 
ages  for  slicy  run  12.  Cross  range  maxi¬ 
mum  and  mean  smear  are  both  success¬ 
fully  reduced.  This  run  does  not  have 
the  largest  amount  of  smear  reduction, 
but  overall  smear  reduction  is  more  read¬ 
ily  evident  than  for  other  runs  due  to  its 
relatively  low  peak  value. 


Figure  5.38:  The  difference  in  cross 
range  smear  between  the  9-aspect  fil¬ 
tered  and  original  slicy  run  12  ISAR  im¬ 
ages.  This  plot  is  the  difference  between 
the  two  smear  functions  in  Fig.  5.37. 
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Figure  5.39:  The  slicy  ISAR  image  generated  with  the  pa¬ 
rameters  given  for  run  12.  This  image  is  unfiltered. 
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Figure  5.40:  The  slicy  ISAR  image  generated  with  the  pa¬ 
rameters  given  for  run  12.  This  image  is  filtered  with  a  3- aspect 
sliding  window,  decreasing  maximum  cross  range  smear  by  ap¬ 
proximately  45%. 
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Figure  5.41:  The  slicy  ISAR  image  generated  with  the  pa¬ 
rameters  given  for  run  12.  This  image  is  filtered  with  a  9- aspect 
sliding  window,  decreasing  maximum  cross  range  smear  by  ap¬ 
proximately  80%. 
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Figure  5.42:  Cross  range  smear  in  the 
original  and  3-aspect  filtered  ISAR  im¬ 
ages  for  slicy  run  1.  Given  this  geom¬ 
etry,  cross  range  smear  is  just  slightly 
increased,  the  only  run  to  do  so.  The 
amount  of  change  is  very  small,  though, 
seen  by  the  overlapping  plots. 


Figure  5.43:  The  difference  in  cross 
range  smear  between  the  3-aspect  fil¬ 
tered  and  original  slicy  run  1  ISAR  im¬ 
ages.  This  plot  is  the  difference  between 
the  two  smear  functions  in  Fig.  5.42. 
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Figure  5.44:  Cross  range  smear  in  the 
original  and  9-aspect  filtered  ISAR  im¬ 
ages  for  slicy  run  1.  Given  this  geom¬ 
etry,  cross  range  smear  is  just  slightly 
increased,  the  only  run  to  do  so.  The 
amount  of  change  is  very  small,  though, 
seen  by  the  overlapping  plots. 
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Figure  5.45:  The  difference  in  cross 
range  smear  between  the  9-aspect  fil¬ 
tered  and  original  slicy  run  1  ISAR  im¬ 
ages.  This  plot  is  the  difference  between 
the  two  smear  functions  in  Fig.  5.44. 
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Figure  5.46:  The  slicy  ISAR  image  generated  with  the  pa¬ 
rameters  given  for  run  1.  This  image  is  unfiltered. 
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Figure  5.47:  The  slicy  ISAR  image  generated  with  the  pa¬ 
rameters  given  for  run  1.  This  image  is  filtered  with  a  3-aspect 
sliding  window,  increasing  maximum  cross  range  smear  by  ap¬ 
proximately  0.3%. 
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Figure  5.48:  The  slicy  ISAR  image  generated  with  the  pa¬ 
rameters  given  for  run  1.  This  image  is  filtered  with  a  9-aspect 
sliding  window,  increasing  maximum  cross  range  smear  by  ap¬ 
proximately  2%. 

5.5  Conclusions 

Overall,  the  cross  range  smear  reduction  technique  as  proposed  by  Bhalla  and 
Ling  [2]  performs  well.  However,  if  the  sliding  window  is  expanded  to  9  aspects,  the 
results  are  more  unpredictable. 

Using  the  originally  proposed  3-aspect  sliding  window,  smear  is  reduced  for 
all  but  a  very  few  cases,  regardless  of  geometry.  In  most  cases,  the  amount  of 
smear  reduction  is  significant,  both  visually  and  by  numerical  analysis.  However, 
the  amount  of  smear  reduction  varies  by  geometry  and,  to  a  lesser  extent,  other 
Xpatch®  parameter  values.  The  smear  reduction  algorithm  also  reduces  the  ISAR 
image’s  overall  scattering  response  level. 

If  the  sliding  window  is  expanded  to  9  aspects,  smear  is  further  reduced  in 
almost  all  cases.  However,  this  smear  reduction  also  reduces  the  overall  response 
level.  Further  increases  in  window  size  would  further  reduce  the  image  brightness. 
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As  can  be  seen  in  the  ISAR  images  that  show  good  smear  reduction,  especially 
for  the  ridged  plate  and  slicy  geometries,  this  smear  reduction  technique  can  per¬ 
form  very  well.  In  most  3-aspect  filtered  cases  it  performs  well,  with  little  image 
degradation.  The  9-aspect  window  usually  produces  even  better  results,  at  a  greater 
reduction  in  brightness. 

Unfortunately,  the  cross  range  smear  reduction  technique’s  performance  is  not 
a  linear  function  of  the  amount  of  cross  range  smear.  Cross  range  smear  is  reduced  by 
different  amounts  in  every  run  and  geometry  combination,  making  smear  reduction 
prediction  very  difficult. 


VI.  Conclusions 


6. 1  Introduction 

This  section  presents  general  conclusions  about  cross  range  smear  in  Xpatch® 
Inverse  Synthetic  Aperture  Radar  (ISAR)  images.  Specifically,  the  two  primary 
facets  of  this  research  are  discussed:  cross  range  smear  characterization  and  reduc¬ 
tion.  In  addition,  future  work  to  characterize  and  minimize  cross  range  smear  is 
suggested. 

6.2  General  Discussion 

Cross  range  smear  is  a  non-physical,  and  therefore  erroneous,  effect  that  occurs 
in  Xpatch®  ISAR  images  generated  with  frequency  domain  computations.  The 
purpose  of  this  work  is  to  facilitate  a  better  understanding  of  cross  range  smear, 
in  terms  of  causes  and  mitigation.  This  work  characterizes  cross  range  smear  in 
terms  of  Xpatch® ’s  input  parameters  and  then  analyzes  a  proposed  smear  reduction 
technique. 

Even  though  cross  range  smear  is  primarily  a  function  of  scattering  geometry, 
parameter  values  have  an  additional  and  measurable  effect.  Through  the  analy¬ 
sis  presented  in  Chap.  IV,  the  following  observations  appear  consistent  through  all 
analyzed  geometries. 

•  Frequency  center  and  bandwidth  have  a  significant  first-order  effect  on  cross 
range  smear. 

•  Ray  density  does  not  have  a  significant  first-order,  linear  effect  on  cross  range 
smear.  Ray  density  has  a  quadratic  impact  over  the  range  of  values  analyzed 
in  this  work.  It  also  has  significant  second-order  influence,  especially  with  the 
angular  increment  parameter. 
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•  The  choice  of  Xpatch®’s  first  bounce  algorithm  or  diffraction  edge  hie  has 

negligible  effect  on  cross  range  smear. 

•  The  scattered  held  computations  took  a  very  long  time,  even  on  a  512-processor 

supercomputer. 

As  seen  in  Chap.  V,  the  smear  reduction  technique  proposed  by  Bhalla  and 
Ling  [2]  performs  well  over  the  geometries  analyzed.  Using  their  3-aspect  sliding  win¬ 
dow,  smear  is  successfully  reduced  in  all  but  a  few  cases.  Often  the  smear  reduction 
is  significant,  with  visually  apparent  improvement.  Unfortunately,  smear  reduction 
also  has  the  effect  of  reducing  peak  response  levels  throughout  the  image,  slightly 
reducing  overall  image  contrast. 

When  the  sliding  window  is  increased  to  9  aspects,  three  times  that  proposed 
by  Bhalla  and  Ling,  cross  range  smear  is  generally  further  reduced.  This  reduction 
comes  with  further  peak  response  level  reduction. 

6. 3  Future  Work 

This  research  represents  a  beginning  of  a  detailed  and  analytic  look  into  the 
cross  range  smear.  As  such,  many  opportunities  for  future  work  exist.  This  work 
concentrates  on  extending  the  smear  characterization  and  testing  smear  reduction 
implementations. 

First,  a  not-yet-released  version  of  Xpatch®  implements  a  cross  range  smear 
reduction  technique  similar  to  that  proposed  by  Bhalla  and  Ling  [2],  When  this 
version  becomes  available,  a  study  of  smear  reduction  similar  to  the  work  performed 
and  presented  in  Chap.  V  is  necessary. 

Second,  more  geometries  with  different  scattering  characteristics  need  to  be 
examined.  Specifically,  cavities  and  occultating  surfaces  present  very  interesting 
targets  with  complex  scattering  properties.  These  targets  have  different  smearing 
characteristics  for  comparison  against  the  targets  analyzed  in  this  research. 
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